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Abstract 

The lepton asymmetry generated by the out-of-equihbrium decays of heavy 
Majorana neutrinos with a quasi-degenerate mass spectrum is resonantly 
enhanced. In this work, we study this scenario within a first-principle 
approach. The quantum field theoretical treatment is applicable for mass 
splittings of the order of the width of the Majorana neutrinos, for which 
the enhancement is maximally large. The non-equilibrium evolution of 
the mixing Majorana neutrino fields is described by a formal analytical 
solution of the Kadanoff-Baym equations, that is obtained by neglecting 
the back-reaction. Based on this solution, we derive approximate analyt- 
ical expressions for the generated asymmetry and compare them to the 
Boltzmann result. We find that the resonant enhancement obtained from 
the Kadanoff-Baym approach is smaller compared to the Boltzmann ap- 
proach, due to additional contributions that describe coherent transitions 
between the Majorana neutrino species. We also discuss corrections to 
the masses and widths of the degenerate pair of Majorana neutrinos that 
are relevant for very small mass splitting, and compare the approximate 
analytical result for the lepton asymmetry with numerical results. 



1 Introduction 



The almost complete absence of antimatter on Earth, in the solar system and 
in hadronic cosmic rays suggests that the universe is baryonically asymmet- 
ric. This conclusion is confirmed by experimental data on the abundances of 
the light elements [I] and precise measurements of the cosmic microwave back- 
ground spectrum [2l |3] . The baryon asymmetry of the universe can be generated 
dynamically provided the three Sakharov conditions [3] are fulfilled in the early 
universe: violation of baryon (or baryon minus lepton) number; violation of C 
and CP] and deviation from thermal equilibrium. Although the Sakharov condi- 
tions are fulfilled in the Standard Model of particle physics (SM) , the smallness 
of the CP- violation in the quark sector and the fact that, given the current 
bounds on the Higgs mass, the electroweak transition is not of first order, do 
not allow for the generation of an asymmetry comparable to the observed one. 
Therefore, the question about the origin of this asymmetry represents a major 
puzzle of modern physics. 

From the theoretical point of view a very attractive explanation of the observed 
baryon asymmetry of the Universe is provided by the baryogenesis via leptoge- 
nesis scenario [5]. In this scenario the Standard Model is supplemented by at 
least two right-handed Majorana neutrinos Ni which couple to leptons and the 
Higgs: 

C = CsM + \N,{i$ ~ M,)N, - KMPrN^ ~ hlN,4>^PLea. , (1) 

where (f> and £ are the Higgs and lepton doublets and <f> = i(J2(t>* ■ We work in 
the mass eigenbasis of the Majorana neutrinos, M = diag(Afi), which are CP 
eigenstates with parities rji = ±1, iVf = tjiNi. The Majorana condition implies 
in particular that N and N are not independent but related by Ni = rjiNfC 
where C = 17270 is the charge-conjugation matrix. This minimal extension of 
the SM can explain naturally the observed small mass scale of the active neutri- 
nos via the famous see-saw mechanism. Furthermore, CP-violating decays are 
responsible for the generation of a lepton asymmetry, i.e. for leptogencsis, which 
is converted into the observed baryon asymmetry by the sphalcron processes [5] . 
The CP-violating parameter receives contributions from the vertex [5] and self- 
energy [2l HI [9l [To] diagrams. If the masses of the heavy neutrinos are strongly 
hierarchical, then the two contributions are comparable. On the other hand, 
if the mass spectrum is quasidegenerate then the self-energy contribution is 
resonantly enhanced and becomes considerably larger than the one from the 
vertex diagram. It has been argued that the resonant enhancement of the CP- 
violating parameter could be compatible with a scale of leptogencsis of the 
order of '~ 1 TeV [lOl |TT] . The lightness of the right-handed neutrino masses 
can in principle induce a detectable non-unitarity of the mixing matrix of active 
neutrinos |12U13j which is very interesting from the theoretical and experimental 
point of view. In addition, a model with right-handed neutrinos lighter than 
10^ GeV does not suffer from the hierarchy problem [T3]. 

The self-energy contribution to the CP- violating parameter has been extensively 
investigated in the literature. One of the first estimates was made in [15j using 
the effective Hamiltonian approach known from kaon physics. For a quaside- 
generate mass spectrum of the heavy neutrinos, M2 — Mi, the obtained result 
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for the CP-violating parameters can be represented in the form: 
_ lm[{h^h)l] KhM,{Mf-M]) 

where i,j ~ 1,2 are flavor indices, and the resonant enhancement is described by 
the function R. For A = 0, the function R would diverge in the degenerate limit. 
Therefore A can be considered as a "regulator", that arises from the fact that 
the Majorana neutrinos are not strictly on-shcll because of their finite lifetime. 
In Ref. [T5] it was found 

^ = 3h^<h^h)MMi + A'hf . (3) 

The regulator is relevant for determining the maximal possible resonant en- 
hancement, which occurs for M| — Ml = ±A, and is given by 

7? - 

J^max — 2\A\ ' 

Note that ^ is proportional to the product M1M2 — A/j ) lm{h^h)l2 which, 
similarly to the Jarlskog determinant, is an invariant (under field rotations) 
quantity characterizing CP-violation in the leptonic sector. In particular, if the 
heavy neutrino masses are equal, i.e. if M'} = M|, then the CP- violating phases 
in the Yukawa couplings can be rotated away and the CP-violating parameter 
vanishes. Extending the formalism developed in (TB] the authors of [T7| studied 
the effects of particle mixing for scalar fields with arbitrary mass splitting and 
also found a considerable enhancement of the asymmetries when the masses of 
the mixed states arc comparable. For decreasing mass splitting the maximal 
enhancement is achieved for — ~ F^^ + 1^22: where Tij denotes the 
absorptive part of the one-loop self-energy diagram. 

An approach based on an expansion of full resummcd propagators around their 
poles was developed and applied to the analysis of Icptogenesis in the SM in 
[21 [HI HH Hn] ■ A related approach based on the diagonalization of the resummed 
propagator was developed in Ref. [50] . The starting point of both approaches 
is the solution of the Schwinger-Dyson equation for the full propagator of the 
mixing species in vacuum. Then the resummed propagators are substituted into 
amplitudes of the lepton-number violating scattering processes, for which only 
stable particles appear in the asymptotic in and out states. The arising expres- 
sions are then interpreted in terms of effective decay widths of the Majorana 
neutrino species and CP-violating parameters are extracted. Within the former 
approach, the decay rates are obtained by extracting the on-shcll contributions 
of the Majorana neutrino propagator to the scattering processes. The resulting 
CP- violating parameter has the form of Eq. ([2]) but with a "regulator" given by 

A = MiTj . (5) 

The result obtained for the CP- violating parameter based on the diagonalization 
approach [5D] also has the form of Eq. ^ , but with 

A = M,T, - MjTj . (6) 

As has been noted in Ref. [H], the difference can be traced back to the way 
how the CP-violating decay rates are extracted from the scattering amplitudes. 
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Furthermore, it has been argued that both the pole mass expansion and diago- 
nahzation approach can be reconciled, and yield a regulator of the form ^ . In 
general there exist two potential sources that can contribute to the regulator, 
arising from the finite width of the two Majorana neutrino species required for 
a CP-violating contribution to the decay rate. Therefore it seems plausible that 
the regulator depends on the widths of both species. As is evident from Eq. 
the latter result predicts a larger value of the CP-violating parameter, especially 
if Ti ^ Tj. A comparison of the various approaches to the calculation of the 
CP-violating parameter in the resonant regime can also be found in |22| . 

In all of the above works the CP-violating parameter is evaluated in vacuum 
and then used to calculate the asymmetry generated by the decays of the heavy 
neutrinos. In the hot and dense plasma thermal corrections to the decay and 
scattering rates can play an important role. An analysis of Icptogenesis within 
thermal quantum field theory was performed in [231 124j . The influence of ther- 
mal corrections on the dispersion relations was studied in [25j [26l [27]. The 
production rate of right-handed neutrinos in a thermal plasma of SM particles 
has been studied in [551 [H] and in [301 131] for temperatures that are large or 
small compared to their mass, respectively. 

The resonant enhancement of the CP-violating parameter becomes maximal 
when the mass difference of the neutrinos is of the order of their width. Ex- 
pressed in a quantum field theoretical language, this means that the spectral 
functions of the neutrinos cannot be described by two distinct peaks, but rather 
by two overlapping resonances. In addition, this implies that it is not sufficient 
to describe the deviation from equilibrium in terms of semi-classical distribu- 
tion functions fNiit) for the two neutrinos i — 1,2. Instead, it is important 
to take a full two-by-two matrix into account, that can describe coherent tran- 
sitions between the neutrino species in terms of cross-correlations. Within a 
purely quantum field theoretical treatment, a natural object to consider are the 
two-point functions {Ni{x)Nj{y)), which possess a matrix structure in flavor 
space. Their time-evolution can be described by the so-called closed-time-path 
(CTP) or Schwinger-Keldysh formalism. Concretely, self-consistent equations of 
motion for the non-equilibrium time-evolution can be obtained by formulating 
the Schwinger-Dyson equation in real time on the CTP. The resulting equa- 
tions are known as Kadanoff-Baym equations. One of the first steps towards 
the analysis of leptogenesis within the formalism of non-equilibrium QFT has 
been performed in [31] where the authors solved the system of Kadanoff-Baym 
equations in the Minkowski space-time using a perturbative expansion. In |33j 
the approach to equilibrium for a scalar field which is coupled to a large ther- 
mal bath was studied. A quantum mechanical calculation of the asymmetry 
generation for a hierarchical mass spectrum of the heavy neutrinos based on 
Kadanoff-Baym equations was presented by the same group in [M] [3S] . We also 
refer to Refs. [36] |37] for previous work on a quantum treatment of leptogene- 
sis. A somewhat different approach was pursued in [35] [Ml UHl HI] as well as 
in [m 1331 133] . In these works the initial system of Kadanoff-Baym equations 
was used to derive Boltzmann-like equations for the quasiparticles and calculate 
their effective in-medium properties like masses, decay widths and CP- violating 
parameters. The vertex and self-energy contributions to the CP-violating pa- 
rameter in a toy model with two real scalar fields and one complex scalar field 
were analyzed in [3S] and [31] . A similar analysis for leptogenesis in the Standard 
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Model was performed in [?5J 23] . Flavor oscillations in the leptonic sector and 
their impact on the generation of an asymmetry were analyzed within nonequi- 
librium quantum field theory in |43] . Oscillations of the right-handed neutrinos 
play a crucial role in the scenario of leptogenesis via oscillations 03 SHI HZ] ■ 
The latter is realized, in particular, in the Neutrino Minimal Standard Model 
(j/MSM) [Tl]. The calculation of the asymmetry generated within this scenario 
is based on so-called density matrix equations [3H]- Another method using the 
perturbative short-time solution of the von Neumann equation and intended to 
check the consistency of the density matrix equations was presented in |49j . 

Summarizing the above, the analysis of leptogenesis in the maximally resonant 
regime poses two problems, which, despite the recent progress, are not yet fully 
resolved. The first one is the enhancement of the CP-violating parameter. The 
existing vacuum calculations lead to rather different results. Furthermore, they 
neglect the medium effects as well as modification of the quasiparticle spectra in 
the hot and dense plasma of the early Universe. It is also a priori unclear whether 
the very notion of quasiparticles is applicable in the maximal resonant regime. 
The second problem is the oscillation of the heavy neutrinos which becomes 
important when the mass difference is sufficiently small. This effect is completely 
neglected in the standard scenario of resonant thermal leptogenesis which, in 
the maximal resonant regime, "merges" with the scenario of leptogenesis via 
oscillations. 

In this work we study the generation of a Icpton asymmetry from first principles 
using an explicit solution of the Kadanoff-Baym equations for mixing Majorana 
neutrinos with an arbitrary mass difference. Following Refs. |34[ I35 j . we will 
treat lepton and Higgs fields as a thermal bath, neglecting the backreaction. 
This simplification has the advantage that it is possible to obtain analytical 
results that can be directly compared to the Boltzmann approach. 
In section [21 we introduce the framework of our computation, explain our no- 
tations, define the quantities that are necessary for a quantum-mechanical ap- 
proach to leptogenesis, and derive the relevant Kadanoff-Baym equations. In 
section [3] we derive an equation of motion for the lepton asymmetry, and discuss 
the generation of an asymmetry from the interaction of the Majorana neutrinos 
with a thermal bath. An analytical solution for the Kadanoff-Baym equations 
and the lepton asymmetry is derived in section 3] using a Breit-Wigner approxi- 
mation. In section [5] we present our results obtained in the Kadanoff-Baym ap- 
proach and compare it to the corresponding Boltzmann approach. In particular, 
we discuss the differences that arise with respect to the resonant enhancement 
of the lepton asymmetry. We also compare the solutions of the Kadanoff-Baym 
equations for various approximations. Finally, we conclude in sectionlHl Details 
of the derivation of the analytical solution can be found in the Appendix. 

2 Closed-time-path formalism 

Leptogenesis is a non-equilibrium process. Within quantum field theory, such 
processes can be conveniently described based on the closed-time-path or Schwinger- 
Keldysh formalism. In particular, we are interested in the non-equilibrium time- 
evolution of the expectation value of physical observables, like the baryon or 
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lepton asymmetry. Consider for example the lepton current operator, 



J£(a:)=^4(a;)7^^a(x). (7) 

a 

We are interested in the time-evolution of its expectation value with respect to 
the physical state of the system, 

= = (^j2^o.ix)rux)^ ■ (8) 

The lepton asymmetry per volume V is then given by 

nL{t) = ^ j^d^xjl{t,x). (9) 

It is the main purpose of this work to compute the time-evolution of the lep- 
ton asymmetry that is generated by the out of equilibrium decays of quasi- 
degenerate Majorana neutrino species within the closed time path formalism. 
In the following, wc will introduce the quantities and equations necessary for 
this computation. 

2.1 Propagators 

Within quantum field theory a system out of equilibrium can be characterized 
completely by its full set of A^-point correlation functions. Since we are inter- 
ested in a temperature regime for which the clectrowcak symmetry is restored, 
the lowest-order correlations are the two-point functions 

^^,,{x,y) = {Tcr{xW\y)). 

Stiix.y) = (TcCW4(y))> (10) 
S'^{x,y) = {TcN,{x)N,{y)) . 

Here a,6 = 1,2 are SU{2) indices, and a, /3 = 1,2,3 as well as i,j = 1,2,3 are 
flavor indices of the leptons and Majorana neutrinos respectively. We use matrix 
notation for the Lorentz-spinor indices. Tc denotes time-ordering with respect to 
the closed time contour, and features the usual minus sign when interchanging 
fermionic operators. In an iS'J7(2)-symmetric state, which we consider in the 
following, one can write 

Se"J{x,y) = S,,S,'^P{x,y). (11) 

Note that the lepton and Majorana neutrino propagators are matrices in flavor 
space, i.e. they encode cross-correlations which are important for coherence 
effects. Throughout this work we assume that the third Majorana neutrino 
is very heavy, and therefore consider the two-by-two matrix in fiavor space 
corresponding to iVi.2. We will frequently use matrix notation for the flavor 
indices. If the two-point functions are known, the lepton current can be easily 
calculated, 

jiix) = -Y,tT[j,S,"'^{x,x)] . (12) 
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Thus, the main task will be to obtain non-equilibrium evolution equations for 
the two-point functions, so-called Kadanoff-Baym equations, and to extract the 
relevant processes which describe the generation of the lepton asymmetry. 
By denoting the time arguments on the chronological branch of the closed-time 
contour by and the one on the anti-chronological branch by x~ , one can 
obtain four propagators with usual real time arguments, 

D++{x,y)^D{x+,V+), D+^{x,y)^D{x+,y-), 
D_+{x,y)=Dix-,y+), D__{x,y) = D{x~ ,y-) , ^ ' 

where D(x,y) stands for any two-point function. The four functions are often 
combined into a two-by-two matrix structure (in addition to the flavor and Dirac 
indices). However, it is easy to see that only two of them are independent. In 
this work, wc prefer to use the more compact notation of jSOj . and decompose 
the two-point function into spectral function Dp and statistical propagator Dp, 

D{x, y) = Dp{x, y) - '-sgndx'' - y'')Dp{x, y) . (14) 

The signum function is either 4-1 or —1 depending on whether x^ or y'^ occur 
'later' on the closed time contour C. From this parameterization it is evident 
that in fact only two independent two-point functions exist. All other two- 
point functions can be expressed as linear combinations. For example, one can 
introduce the Wightman functions D^{x,y) = Dp{x,y) =F ^Dp{x,y). In terms 

of the matrix notation, these would correspond to £)> = D y and £)< — D^ 

Finally, it is useful to introduce retarded and advanced functions, 

DRix,y) = e{x" -y")Dp{x,y), (15) 
DAix,y) = -Qiy"" -x'')Dp{x,y). (16) 

A useful relation for the Majorana neutrino propagator, which results from the 
relation = rjiNi, is 

S'^{x,y)^7nrj,CS'\y,xfC-' , (17) 

implying that 

^F,R/A(^^y) = mV3CS^pj^jj^{y,xfC~'^, 
S^^{x,y) = -rj,r^jCSp{y,xfC-'. 

CP-conjugate propagators 

For leptogenesis, it is useful to consider separately CP-odd and CP-even effects. 
In particular, the generation of the lepton asymmetry is a CP-odd effect. In 
order to discriminate between CP-odd and CP-even, one can perform a CP 
transformation of the system, and determine the transformation properties of 
the two-point functions. 

The Higgs field transforms like (j>{x) — > (l)*{x) under CP, where a; = (a;°,— x). 
Thus the propagator gets transformed in the following way: 

^0ab(^'2^) = (TcM^Wbiy)) ^ {Tc(t)l{x)(j)b{y)) = A^^^{y,x) . 
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Motivated by this observation, we define the CP-conjugate propagator as 



A^^^{x,y)=A^^^{y,x). (18) 
From this, one can deduce the various components introduced above, 

It is important to note that for a system in a CP-symmetric state the propagator 
and the CP-conjugate propagator are equal, A^^^{x,y) = A^^j^{x,y). 

The lepton field transforms like i{x) CPF{x) and i{x) ^ £^ {x)PC ^ 
—f"{x){CP)^^, where P — 70. Accordingly, the propagator gets transformed 
like 

sCabi^^y) {cp),^s/li,syrx)[{cpr%. , 

where we have displayed Dirac indices for clarity. The CP-conjugate propagator 
is then given by 

SC^{x,y) ^ CPSdy,xfiCP)-^ . (19) 

Here we have again used a matrix-notation for the Lorentz-spinor indices, and 
the transpose on the right-hand side refers to these indices. For the components, 
one finds similarly as above 

Slab^.FM/Ai^^y) = '^PS4a$,KA/B.iy'^)^('^Py^ ' 

S"Jp{x,y) = -CPSi:^{y,xf{CP)-K 

For a system in a CP-symmetric state, one has Si"l^{x,y) = Si^^{x,y). 
The CP-conjugate propagator of the Majorana neutrinos is given by 

S'\x,y)=CPS^\y,xf{CP)-^ . (20) 
For the components, the relations are also analogous to the leptons. 



2.2 Kadanoff-Baym equations 

The non-equilibrium equation of motion for the two-point functions can be de- 
rived by a variational principle from the so-called two-particle-irreducible (2PI) 
effective action T[A^, S^, S]. Technically, it is a Legendre transform with re- 
spect to the generating functional in the presence of bi-local sources. All the 
information about the quantum system is encoded in the effective action. 
The equation of motion for the two-point functions follows from the stationarity 
conditions. 

In general, the 2PI effective action can be parameterized in the form 

r[A^,5£,5] =5 + ri^[A^,5f,5]+r2[A^,5f,5], (22) 

where the classical action S = d^xC is the tree-level contribution, T^^ the 
one-loop contribution, and T2 contains diagrams with at least two loops. More 
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Figure 1: Two- loop diagram contributing to the 2PI effective action. 




Figure 2: One-loop diagrams contributing to the lepton (left) and Majorana 
neutrino (middle and right) self-energies. 



precisely, r2[A0,5^,5] is equal to the sum of all 2PI ('skeleton') Feynman di- 
agrams which have no external legs. Furthermore, the internal lines represent 
the full propagators A^jS^, and S. The two- loop contributions are shown in 
figure [TJ 

The resulting equations have the form of Schwinger- Dyson equations, evaluated 
on the closed-time path. For the leptons and the Majorana neutrinos, they read 
(suppressing flavor, SU(2) and Lorcntz-spinor indices for the moment) 

|I=0 O Si-\x,y) = Sio\x,y)-J:i{x,y) , (23a) 

^=0 ^ S-\x,y) = So\x,y)~^Nix,y) . (23b) 

Here Se^^{x,y) is the inverse of the full lepton two-point function (i.e. the 
inverse of the full propagator), SiQ^{x, y) is the inverse of the free lepton prop- 
agatoi0, 

iSia\x,y) = t0JabS"^6c{x - y)PL , (24) 

and (x, y) is the self-energy given by the functional derivative of the effective 
action r2[A0, Sg,S], 

^tt.abi^^ v) = pl^] — : ■ (25) 

Graphically, this corresponds to opening one lepton line of the diagrams con- 
tributing to Thus, for example the two- loop contribution to r2 yields 
the one-loop self-energy. The one-loop contribution to the lepton self-energy 
is shown in figure [5] We stress again that the internal lines of the diagrams 
represent the full propagators. This means that the equation of motion is a 
self-consistent Schwinger-Dyson equation. 



^Strictly speaking, in Dirac notation the inverse exists only when restricting the spinors to 
the left-handed sub-space. Nevertheless it is possible to work in Dirac notation when properly 
inserting left-handed projectors , so that effectively only the left-handed parts contribute. 
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Similarly, the free inverse propagator for the Majorana neutrinos reads 

iSo\x, y) = (z4 - M,) 5'^c{x - y) , (26) 

Later on, it will be necessary to cancel divergent contributions to the self-energy 
of the Majorana neutrinos. Therefore, we introduce suitable mass and field 
counterterms. Their most general form can be obtained by applying a field 
rescalingiVi {ZijPL+r]iZ*jr]jPfi)Nj compatible with the Majorana conditions 
and using the most general form of the mass matrix for the mass counterterms, 
5Mij = i]iSMjirij. The resulting counterterms are given by pUl \5T\ 

6C = ^Nl^^{dZ^ + SZ)PB,0+^{'nSZ'^r] + 'n5Z*ij)PL0-^{5Z^M + 2dM^ 

+ MriSZ*T])PR - i {rjSZ'^TjM + 25M + M5Z)Pl^N , (27) 

The (renormalized) self-energy has two contributions, 

^N%{x,y) ^ + 5^N%{x,y) ■ (28) 

The latter arises from the counterterms (P7|) . 

5Y.n'^ (x, y) = -i5^n'' {x)5c [x - y) , (29) 

~[^{SZ^ + 5Z)Pr$,^ + ^{7]5Z'^r) + r^5Z*T^)PL$x - i(,5ZtM 
+ 2(5 Aft + MT]SZ*ri)PR - KtjSZ'^tjM + 26M + MSZ)Pl^ , (30) 

The Majorana self-energy inherits the relation 

j:N'\x,y) = 7j,7j,CJ:N'\y,xfc-^ , (3i) 

of the Majorana propagator. Note that also the counterterm-part fulfills this 
relation. 

The equation of motion can be brought in a more explicit form by convolut- 
ing it with the full propagator, and inserting the decomposition into statistical 
and spectral components. Then one obtains a system of two coupled integro- 
differential equations, so-called Kadanoff-Baym equations. 







z45,f(x,y) = [ d^zJ:e;''{x,z)Sefiz,y) 

d^zi:,'P{x,z)S,f{z,y), (32a) 
t0,S,f{x,y) = I d^zJ:e''^'>{x,z)S,f{z,y). (32b) 

JyO 

The corresponding Kadanoff-Baym equation for the Majorana neutrinos reads 
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(in matrix notation for the flavor indices) 

- M - 5Y.N{x))SF{x,y) = / dt^ zY.^ p{x, z)SF{z,y) 

JQ 

d'^zJ:NF{x,z)Sp{z,y) , (33a) 
{i0.^ - M - ST,N{x))Sp{x,y) = / d^zJ:Np{x,z)Sp{z,y) . (33b) 

These equations contain 'memory integrals' on the right-hand side, which inte- 
grate over the history of the system. Thus, the change of the two-point func- 
tion at a given instant of time depends on the entire history of the system, 
i.e. the equations are non-Markovian. Furthermore, the equations explicitly 
depend on the two time-arguments and j/" characterizing the two-point func- 
tion. For a system in equilibrium, they would only depend on the difference 
gO — 2;° — y". Out of equilibrium, the system also depends on the average time 
X° — (a;° + ?/")/2. For this reason, the equations in the form given above are 
sometimes referred to as two-time Kadanoff-Baym equations. We stress that 
these equations are valid also far from equilibrium. For systems which are close 
to equilibrium, one expects that the two-point functions depend only weakly 
on X". In this case, it is customary to perform a gradient expansion, which 
yields Markovian quantum kinetic equations. In addition, one may perform a 
Fourier transformation with respect to (Wigner transformation). Then, a 
very common further simplification consists of making a particular ansatz for 
the frequency-dependence of the two-point functions, including a Breit- Wigner 
ansatz or a quasi-particle ansatz. The latter is known to lead to Boltzmann-type 
equations. 

On the other hand, it is also possible to follow a different strategy, namely to 
retain the non-Markovian, two-time structure of the equations. If, in addition, 
the self-energies are computed using the complete propagators, the resulting 
equations are non-linear integro-differential equations. Equations of this type 
have been solved numerically with increasing complexity over the past years, 
and it has been shown that they can describe the process of kinetic and chemi- 
cal quantum thermalization for systems that are initially far from equilibrium, 
see e.g. [SH [521 EH ESI EEl EZl EH]- A significant simplification can be achieved 
when treating the self-energies as time-translation invariant functions, which 
reduces the Kadanoff-Baym equations to linear integro-differential equations. 
Physically, this can occur when the quantum field that is out of equilibrium is 
weakly interacting with a thermal bath, and the interaction can be described 
by self-energies that solely depend on propagators associated to particle species 
belonging to the thermal bath. The approximation of time-translation invari- 
ant self-energies then corresponds to neglecting the backreaction of the non- 
equilibrium quantum field on the thermal bath (more precisely, to neglect the 
possible deviation from equilibrium of the self-energies that is caused by the 
backreaction). For that case, approximate analytical solutions of the Kadanoff- 
Baym equations within the two-time representation have been derived in |33j . 
Moreover, these solutions turn out to be well-suited for describing the evolution 
of the Majorana neutrino propagator, and have been used to study Icptogenesis 
in the hierarchical limit |34[|35j . The reason is twofold: First, the corresponding 
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self-energies depend on the lepton and Higgs propagators, and the latter belong 
to the thermal bath of 0(100) SM particles that are kept in equilibrium due to 
gauge interactions. Second, the Majorana neutrino itself interacts weakly with 
the SM via the Yukawa couplings. In the present work, we will follow the same 
strategy to obtain analytical solutions of the Kadanoff-Baym equations for the 
resonant case. 

Although the two-time Kadanoff-Baym equations shown above are valid far 
from equilibrium and do not require any a priori assumptions of the spectral 
properties, they do have a restriction which is not so commonly emphasized. 
Namely, the equations in the above form implicitly assume that the initial state 
of the system, corresponding to a;" = j/*^ = 0, is a Gaussian initial state. This 
means that the connected part of all higher correlation functions (three-, four- 
point correlations and so on) vanish at the initial time. Physically, this means 
that the interactions are 'switched on' just at the time t = 0. However, it 
is possible to generalize the Kadanoff-Baym equations such that they can also 
describe systems with more realistic initial conditions, if necessary [59]. Below, 
we will in fact include higher correlations in an effective way by taking the 
evolution of the system before the 'initial' time into account. 

2.3 Self-energies 

The self-energies, which enter into the Kadanoff-Baym equations on the right- 
hand side, contain all the information about the interactions of the system. As 
discussed above, the one-loop self-energy is obtained from the two-loop contri- 
bution to the 2PI effective action (see figure [T|), 

Tl^[A^,S,,S] - y^hlhfi, f d^x f d''ytr\s^\y,x)PLS,:^{x,y)PR 

xe''^e^^A^Jx,y). (34) 

The one-loop lepton self-energy obtained from differentiating Fl'^ with respect 
to the full lepton propagator Si{y,x) reads (see first diagram in figure 

^Cb(^^y)iL = h^,h]pPRS'^{x,y)PLe-'e-''A^^^{y,x) , (35) 

where summation over i,j is implied. For a SU(2)-symmetric state, one can 
insert that A^^j^ = SabA^ and S^"^ = SabSi"^. Then the self-energy is itself 
proportional to Sab, and one can write S^"^ = Sab'^i"'^ , with 

E,"'3(a;,y)iL = -h^,hlpPRS'^x,y)PLA4y,x) . (36) 

The one-loop contribution to the self-energy of the Majorana neutrinos reads 

^N'{x,y)iL = hlh^,PLSrJ{x,y)PRe''\^^A^^^{x,y) 

r^.V,hlh},^PPLSCc{S:,y)PRP^'"'e^'^Hdi^,y) ■ (37) 

Here summation over a, /3 is implied. After contraction of SU(2) indices 

Ejv^^(a;,y)iL = -2hlhpjPLSe^"^{x,y)PR 

- 2ij,ij,hlh;^PPLSe4,"^{x, y)PRP , (38) 
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where we have introduced an abbreviation for the lepton-Higgs loop 

S,^''\x,y) ^ CPSi/"iy,xfiCP)-' ^S"\x,y)A^{x,y) . (39) 
Similar to the propagators, it is useful to define CP-conjugate self-energies, 
S,"^(x,2/) ^ C7PI]A(y,x)^(CP)-i , 

f:N \x,y) = CPi:N'\y,xf{CP)-' . 

Then one finds for the CP-conjugated one-loop lepton self-energy, 



Sf"''(x,y)iL = -hl^hjpPRS'\x,y)PLA^{y,x). 
and for the Majorana neutrino 

tN'\x,y)iL = -2hJ'^hp^PLSi^"^{x,y)PR 

- 2r],r],hlhp,PPLSi4,''^ix, y)PRP . 



(40) 
(41) 



(42) 



3 Equation of motion for the lepton asymmetry 

The most important observable for leptogenesis is the lepton asymmetry. Its 
time-evolution is usually computed by subtracting the lepton and anti-lepton 
particle densities. However, since particle number distributions are a classical 
concept, we have to replace them by a suitable quantum-mechanical observable. 
A suitable starting point is the operator 



(43) 



which corresponds to the Noether-current for the lepton fiavor a for a ~ l3. We 
further define a time-dependent lepton-number matrix by 



(44) 



The entries on the diagonal correspond to the lepton asymmetry in fiavor a at 
time t per volume V. The total lepton asymmetry per volume summed over all 
fiavors is 

a 

which is independent of the choice of the flavor basis. The expectation value of 
the lepton current operator can be related to the lepton propagator. 



jLpi^^) ^ (JL^i^)) = (4(a:)7%(x)) = -tr 



(46) 



In the presence of lepton-number violating processes, its divergence can be non- 
zero. In terms of the propagator, it can be written as 



= -tr 
= -tr 



j,id^,+di')Se"('{x,y) 



r,Se''^{x,y) + Se'^^{x,y) 



x=y 



= -J rf*2tr 



{x, z)Se''^{z, y) - Sr{x, z)S,^^'\z, y) (47) 
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-tr 



= -tr 



(.x,y)$y 



V ' 



x=y 



= - / d^ztr 



Si 



{x, z)Sf^{z, y) - Se^'^ix, z)SeQ ^^''(z, y) 



x=y 



(48) 



In the last line we have inserted the free inverse lepton propagator ^" (x, y) = 
6°'^0x^cix — y)PL and used the chiral symmetry relations Sf""^ = PlSi"^ = 
St"'^ Pr- Now, one can replace the free inverse propagator Si^^ on the right- 
hand side by using the Schwinger-Dyson equation Sg'^ = Si^^ — . The parts 
containing Sg~^Sg do cancel, and we obtain the following equation of motion 
for the lepton current: 

^^^fLc.B{^) = - I d^ziAi:r{x,z)SpP{z,x)-Sr[x,z)i:P^{z,x)\ . (49) 



From this, one can immediately derive an equation of motion for the lepton 
asymmetry (assuming a spatially homogeneous system for which d^j^ = doj^), 



dt 



d X I d yti- 



^r{x, y)Spf'{y, x) - Sr{x, y)I],^^(y, x) 



(50) 

where x = (t,x). This can be considered as the master equation for the 
quantum-mechanical treatment of leptogenesis. For future reference, we also 
remark that one can rewrite the equation equivalently as 

driLapit) 



dt 



-i / d^x 



d'^yti- 



i:t''{x,y)Spf'{y,x) - i:r{x,y)Sr{y,x) 



(51) 

where we have used that tiYi/^ (x,y)Si^°' (jj,x) = tiSg"'^ {x,y)Yip^ {y,x), and 
substituted x — x and y — !• — y in the spatial integrals for the second term. 
The 'master' equation ()50p can be brought in a more explicit form by inserting 
the decomposition into statistical and spectral components. In addition, for 
a spatially homogeneous system the two-point functions depend only on the 
difference of the spatial coordinates x — y. Then it is convenient to switch to 
Fourier space in the spatial coordinates. 



dt 



dt' 



X tr 



d^p 
(2^ 

i:,'^/^{tX)Sif^it',t) - Ee'p^{t,t')Se;l{t',t) 



5,«^(t,t')S.?^';(i',i) + 5.7p(i,^')S.^^(i',0 



(52) 



By integrating the upper equation with respect to t, one can easily obtain the 
time-evolution of the lepton asymmetry riL^pit), 



X tr 



dt' / dt" 



d^p 
(2^ 



j:r,l{t',t")Sef^{t",t') - ^tF\{t\t")S,f^{t" ,t') 



s,';i{t',t")^,f^{t".t') + s,f'^{t',t")^i;i{t",t') 



(53) 
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If the lepton propagator is approximately flavor-diagonal, S'^"^ = S^^Si, it is 
easy to see that the integrand is symmetric with respect to t' ^ t" . Using that 

/o dt' /o' dt"{f{t', t") + f{t", t')) = /q dt' /o dt"f{t', t"), it is possible to obtain 
a simplified equation, 

^ (27r 



nLc.p{t)=^ I 7^ fdt' f dt"ti-\j:efp{t'X)SeFp{t",t') 
' ] Jo '- 

-i:if^{t\t")S,ppit",t')\ . (54) 



This equation essentially agrees with Eq. (20) of However, we find it more 
convenient to use Eq. ((5T|) in the following. 



3.1 Generation of an asymmetry 

In the following, we will explicitly compute the time-evolution of the lepton 
asymmetry in the full two-time approach using the two-loop truncation of the 
2PI efl[ective action. For simplicity, we will neglect washout effects and neglect 
the expansion of the universe for the moment. Furthermore, we will assume 
that the generated asymmetry is small, such that one can expand around a CP- 
symmetric state. For that purpose, we parameterize the two-point functions 
and their CP-conjugates as 

Dix,y) = D^x,y) + ^SD{x,y), (55) 

Dix,y) = D'ix,y)~^SD{x,y), (56) 

where D € {A^,^^}. For a CP-symmetric state (i.e. with no lepton asymme- 
try), 6D vanishes. Inserting this parameterization into the equation of motion 
(j5ip for the lepton asymmetry, and expanding in 6D, yields at leading order 

= ^1.J d'^j^d^ir [(l],"^(x,y) - s/"(x,2/)) \sS,^{y,x) 

(57) 

where we have assumed 81°"^ ~ S^^Si, and the propagators entering the self- 
energies are to be evaluated with the symmetric parts of the propagators, I?*, 
as indicated by the subscript. The washout terms would contribute at next-to 
leading order in SD. 

Now we are ready to insert the explicit one-loop expressions for the lepton self- 
energy. This yields the following equation for the lepton asymmetry: 



dt V 



d^x / d yli 



Pr [S''{x,y) - S'\x,y)) PLSi^{y,x) 



(58) 

For brevity, we have dropped the superscript 's' on lepton and Higgs propaga- 
tors, and used the abbreviation Si,f,{x, y) = Si{x, y)/S.(f,{x, y) for the lepton-Higgs 
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loop. By performing manipulations similar to the ones above, 



f (Pp f* /■* 

n-Lapit) = -ihath]^ J dt' i 



df 



tr 



This equation provides a quantum field theoretical expression for the lepton 
asymmetry that is generated from the wave or self-energy type contributions. It 
is appropriate to study the resonant case, where these contributions dominate. 
As input, it requires the two-point functions for the lepton, Higgs and Majorana 
neutrino fields. Since lepton and Higgs interact much more strongly than the 
Majorana neutrino, it is often assumed that they arc in thermal equilibrium. 
Thus, one needs to provide (at least) the equilibrium expressions for the lepton 
and Higgs propagators, as well as the non-equilibrium Majorana propagator. In 
particular, the information on the CP violation is encoded in the off-diagonal 
components of the Majorana propagator. In the following, we will evaluate the 
lepton asymmetry for different choices for these propagators, that correspond 
to different approximation levels, and enable a comparison to the conventional 
Boltzmann result. 



3.2 Asymmetry in a thermal bath 

In the following, wc will derive an expression for the lepton asymmetry assuming 
that lepton and Higgs fields can be treated as a thermal bath, and that the back- 
reaction can be neglected. This setup has been considered in Refs. [33 135] for a 
hierarchical neutrino spectrum. Concretely, this means that we assume that the 
two-point functions for lepton and Higgs fields are time-translation invariant, 

Se{x,y) = S/''{x-y), 

A^{x,y) = A^'\x-y). (60) 

In the following we will explore the consequences of this assumption. The 
most important one is that it is possible to obtain an analytical solution of 
the Kadanoff-Baym equations for the Majorana neutrino fields, which will be 
discussed in the following. 

Non-equilibrium Majorana Propagator 

As a first step, we derive a Kadanoff-Baym equation for the retarded and ad- 
vanced propagators 

SR{x,y) = e{x" -y°)Sp{x,y), 

SA{x,y) = -e{y"-x")Sp{x,y). (61) 

First, we apply the differential operator that appears on the left-hand side of 
the Kadanoff-Baym Eq. p3ap to the first line, 

- M - SJ:Nix))SR{x, y) 
= ijo6{x" - y°)Sp{x,y) + e(x° - y")(i4 ~M~ <5Sw(x))5p(x, y) . 
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Here we have used that the spectral function at equal times is determined by 
the commutation relations of the Majorana neutrino fields, 

S'^{x,v%o=y. = i{{N,{x),N,[v)])x<>=yO = i7o'5*^''5(x - y) . (62) 

Using the Kadanoff-Baym Eq. (|33ap for the spectral function, wc obtain 

- M - SY,n{x))Sr{x, y) 

= ^S''6{x-y) + e{x° -y'') d^zi:^ p{x,z)Sp{z,y) . (63) 

JyO 

Next we rewrite the memory integral by using the identity 

d^z = J d^z [e(x" - z")e{z'> - y") - e(y" - z°)e(z° - x")] . 

where J d'^z = Jd^z denotes a time- integration over the whole real axis 

and over spatial coordinates. Due to the factor Q{x'^ — y^) in Eq. only the 
first term contributes. Inserting this relation yields a Kadanoff-Baym equation 
for the retarded propagator. The equation for the advanced propagator can be 
derived analogously. They can be summarized as 

- M - SY.N{x))Sii(^A){x,y) 

= -S''5{x-y) + j d'^zT.NR{A){x,z)SR(A){z,y) . (64) 

For thermal lepton and Higgs fields, the self-energies are translation invariant 
and depend only on the relative coordinate s = x — y (note that 5T,m{x) = 
5T,m{x + a)). This implies that the equations for retarded and advanced propa- 
gators become invariant under the translation x ^ x' = x -\- a, y ^ y' = y + a. 
Furthermore, their initial conditions are fixed by the equal-time commutation 
relations, as for the spectral function. Consequently, the retarded and advanced 
propagators themselves are translation invariant, 5'i?(A) (2;, y) = (2;', 2/')- 

This means that they can depend only on the relative coordinate as well, i.e. 
Sji(A) (x, y) = SfK^A) {x — y)- For functions that depend on x — y, it is convenient 
to switch to momentum space. We define the Fourier transformation by 

^^""-y^^ J ^"''^'"'^^(p). (65) 

The Fourier transformed equation for retarded and advanced propagators reads 
i - Af,) S^^ S^n'' {p) - J:^%a) (p)] ip) = -S'' ■ (66) 

This is an algebraic equation, with Dirac as well as flavor matrix structure. The 
interaction with the thermal bath of lepton and Higgs is encoded in the self- 
energy. The sum of the retarded or advanced self-energy and the counterterm 
contribution yields renormalized expressions for the self-energy. 

Now we can return to the Kadanoff-Baym equations for the spectral and statis- 
tical Majorana propagators, Eq. p3ap . The spectral function is directly related 
to the retarded and advanced propagators, 

Sp{x - 2/) = Sr{x - y) - Sa{x - y) . (67) 
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Consequently it is also translation invariant. 

The KadanofF-Baym equation for the statistical propagator, Eq. p3ap . can be 
written as 



d 



4, 







- M^) S'^ - sj:n"'{x)) s{x -z)- ^Nt{x,z)] ^■'■(^,2/) 

d^z^N'F{x,z)S''/{z,y) . (68) 



In this form, the homogeneous part is written in the first line and the inhomo- 
geneous part in the second one. The solution can be written as the sum of a 
special solution of the inhomogeneous equations plus the general solution of the 
homogeneous equation. For the former we make the ansatz 

Sp{x,y)inhom = - d'^uS]^{x,u) d'^zY.Np{u,z)S'j{z,y) . (69) 



The application of the operator in square brackets in Eq. ([BS]) will turn the re- 
tarded propagator into a delta function. Then integration over u becomes trivial 
and yields the right-hand side of the Kadanoff-Baym equation, as required. 
Next, we consider the homogeneous equation, i.e. we set the right-hand side 
in Eq. (|68p to zero. The general solution will depend on the initial conditions. 
Since only the first derivative appears, the general solutions can be completely 
specified by the initial value of the propagator itself, Sp (x,y)\,j.o^yO^Q. Here, 
we make the ansatz 

S'^{x,y)ho.n = -Jd\J d\ S]^(x,(Q,u)) A'''{u,v)S'i{{0,^r),y) , (70) 

where ^'^'(u,v) is a free function. Using again the equation of motion for the 
retarded propagator, one can see that this ansatz indeed solves the homogeneous 
equation (except for delta-functions located at the initial time surface; these 
determine the initial conditions but do not influence the validity of the solution) . 
Also, it is important that the dependence on y solves the hermitian conjugated 
Kadanoff-Baym equation, as required. The initial condition to which the ansatz 
corresponds can be extracted by letting x°,y^ from above. In this case we 
can replace the retarded and advanced functions by the spectral function, and 
use its equal-time limit. This yields 

Sp{x,y)hom\x'>=y°=o = «7o^''' (x, y)i7o ■ (71) 

Thus there is a simple one-to-one correspondence between the function A^^ (x, y) 
and the initial condition of the statistical propagator. The general solution is 
given by the sum of both contributions, 

Sj^ (x, y) = S"^ (x, y)inhom + S"^ [x, y)hom ■ (72) 

This is a formal solution, which is valid even when lepton and Higgs were not 
in equilibrium. However, when they are treated as a thermal bath, then the 
retarded and advanced functions as well as the self-energy depend only on the 
time-differences. Thus, in a thermal bath, the integrand in the inhomogeneous 
contribution involves only time-translation invariant functions. Therefore, one 
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might expect that the inhomogeneous contribution is itself time-translation in- 
variant. However, this is not the case because of the lower integration limits, 
which start at the initial time t — 0. The deeper reason behind this is that 
standard Kadanoff-Baym equations rely on a hidden assumption. Namely, it is 
implicitly assumed that all the higher correlation functions (three-point, four- 
point function etc.) vanish at t = 0, i.e. that the system is in a Gaussian state 
at t = 0. Then, higher correlations have to build up, which is reflected in the 
finite integration range starting at t = 0. Physically, we expect that higher 
correlations are present at any time. A simple way how to incorporate these is 
to consider also the evolution of the system at times 'before' the initial time, i.e. 
for t < 0. Technically, this means that all integrations start at t = — oo instead 
of t = 0, 



inriojn 



Sr{x, u)Y,nf {u, v)S''l (u, y) 



(73) 



Then, the inhomogeneous contribution will no longer vanish at x*^ 



y 



0, but 



will be determined by the pre-evolution of the system (i.e. at times ,y'^ < 0). 
On the other hand, the ansatz Eq. ((70|) still represents a valid solution of the 
homogeneous equation for non-zero time arguments. It describes a Gaussian 
deviation from equilibrium occurring at x^ = y'^ = 0, and can be traced back to 
an external source term, that is also present within the usual Gaussian frame- 
work and that originates from the contribution of the density matrix to the 2PI 
effective action [59l |60] . 

Let us now revisit what happens when lepton and Higgs are treated as a thermal 
bath. Then, since retarded and advanced propagators as well as the self-energy 
depend only on the difference of the time arguments, and the integration range 
is now over the whole real axis, the inhomogeneous solution itself will be time- 
translation invariant, i.e. 



Sp (x, y) inho 



POO 


/•oo 


/ d\ 


/ d\ 


/-oo J 


— oo 



s 



S}^{x~u)j:M'}^iu-v)S'iiv^y) 

ij th I \ 

[x-y] . 



(74) 



This means that, in a thermal bath, the inhomogeneous part simply corresponds 
to the thermal equilibrium propagator. The deviation from equilibrium is there- 
fore described by the homogeneous contribution Eq. (j70p . In particular, from 
Eq. ([7T1) it follows that the function A appearing in Eq. ([70]) specifies the devia- 
tion from equilibrium at the 'initial' time x^ = y^ =^ 0, 

AS']^(a;,?/)Uo^,y0^o = {S''^{x,y)- S''p^\x-y))^a=yQ^Q 

= t7oA^^'(x,y)i7o . (75) 

Since we are interested mainly in spatially homogeneous systems, it is convenient 
to switch to spatial momentum space, but to retain the dependence on time. 
In particular the function A fulfills A*^ (x, y) = A*^ (x — y) . Then the formal 
solution for the statistical propagator can be summarized as 



{t,t') 



du 



dv 



S^i^{t,0)A'^S'i„{0,t') 



(76) 
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As discussed above, in a thermal bath the inhomogeneous part of the solution 
yields the thermal propagator, and the homogeneous part describes the deviation 
from equilibrium. The solution for the statistical propagator can thus be written 
as 

-5i,^pWz7oA5^'p(0,0)*7o5!^p(-t') . (77) 

Due to the retarded and advanced functions, the second line is non-zero only 
when t > and t' > 0. This means that it describes the impact of a distortion 
that has occurred at time t = t' = 0. In other words, it describes how the 
system reacts and how the distortion is propagated to positive times t, t' > 0. 
We emphasize that there were no assumptions necessary that would limit the 
size of the 'distortion' described by AS'|,'p(0, 0). This means that the equation 
is valid also for large deviations from equilibrium, as long as the assumption 
of a thermal bath is justified. For example, for 'zero initial abundance', the 
Majorana propagator at t = t' = is the vacuum one, Spp{0,0) = S'f.'™'^(0). 
This means one has to choose AS'^'p(0, 0) = S'f.'™=(0) - S'f,'*''(0). 

Asymmetry for thermal lepton and Higgs 

Starting point for our calculation of the lepton asymmetry is the quantum field 
theoretical expression Eq. ([5^ . We insert the thermal lepton and Higgs propa- 
gators, and write the Majorana propagator in the form 

S'^^M = S'/;\t t') + AS'^^it,t') . (78) 

Furthermore, we extend the lower limits of integration in Eq. (|59p to — oo in 
order to take higher correlations at the 'initial' time t — into account, as 
explained above. Then the contribution to the asymmetry that involves only 
equilibrium propagators vanishes, because the lepton asymmetry is zero in ther- 
mal equilibrium. Next, we insert the solution Eq. ([77)1 of the Kadanoff-Baym 
equation for the statistical propagator, which can be written as 

A5|p(i,t') = -5fp(t)*7oA5^'p(0,0)z7o5^p(-t') (79) 
= e(i)e(t')5;'p(i)»7oA5|,'p(0,0)z7o5;/p(-i') ■ 

It is non-zero only for t > and t' > 0. This means that the lower limits of 
integration in Eq. ([5^ are again zero. Note that it was nevertheless important 
to extend them in the first place, because otherwise the 'equilibrium' contri- 
bution would not have vanished due to unwanted correlation build-up effects. 
Altogether, we arrive at the following expression for the lepton asymmetry in a 
thermal bath: 



riLc^pit) = iho^.h]^ I f dt' f dt" (80) 



JO 



tr 



Pn [AS'^^{t',n AS';^{t',n) PLSe,^^{t" t') 
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Weakly interacting thermal bath 

A case of particular interest arises when the interactions of Icpton and Higgs 
are strong enough to maintain thermal equilibrium, but at the same time weak 
enough to justify the use of free thermal propagators for their two-point func- 
tions. In the following we will explore the consequences of this assumption, and 
we postpone the discussion whether this assumption is justified within the SM. 
The case of free thermal propagators is also interesting because it enables a 
rather direct comparison with the Boltzmann result. 

The free thermal propagators for the Higgs and lepton fields in the time/momentum 
representation are given by 

1 + 2 f'^'^(k) 

^Mt-t') = ^±11 cos {kit -t')), 

^/Mt't') = isin(fc(t-0), (81) 

and 

S,7^it-t') = i^M!i^|zip:cos(fc(t-t'))-*7osin(fc(t-t')) 

-^Ckl^-*') = |^^sin(fc(t-t')) + «7ocos(fc(t-t'))| , (82) 

where k = |k|, f^\k) = 1/(6*^/^ - 1), and f'^k) = 1/(6*^/^ + 1). The chiral 
propagator is obtained by multiplying with the projector from left. Note that 
the correct chiral structure is nevertheless taken into account by the projectors 
Pl/r contained in the self-energies. In addition, similar to Ref. j35j . we take 
the effect of a thermal width into account qualitatively within the the lepton- 
Higgs-loop Si^^ appearing in Eq. ([50)1 by considering damped Breit-Wigner 
propagators 

A0l.%)k(^-O ^ A,-^),(^-^')e-^^l*-*'l/^ 

^Ap)k(^-^') ^ 5/;(p)k(^-i')e-^'l'-*'l/^ (83) 

Evaluating the lepton-Higgs-loop Sgtpp defined in Eq. with the thermal lep- 
ton and Higgs propagators specified above yields 

Si%\^{t -t') ^ j [sr^^{t - i')A0^%(i - 1') + sf^^it - t')A,;v^ - ^')' 

_0_e-^.,\u\/2i^S,1j,+ e'(^+«)" + e-^('=+')" (84) 

The coefficients introduced in the last line are given by 

'S'4p = ^(fc7o±k7)/£0(fc,g), 

54j = ^(fc7o±k7)(/4g) + /,(fc)) , 
where k = p — q, g = |q| and k = |k|, F^^ = + F^, and 

Mk,q)^l + Uiq)~h{k). 
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4 Analytical Breit-Wigner approach 



The equation for the lepton asymmetry, Eq. (j80|) . together with the expression 
(|79p for the non-equiUbrium evolution of the Majorana neutrinos in flavor space, 
and with the Schwinger- Dyson equations for the retarded and advanced 
neutrino propagators can be solved numerically. Before discussing the numerical 
solution, it is however very useful to pursue an analytical approach. This is the 
content of this section. 

We will proceed in three steps. First, we compute the retarded and advanced 
two-point functions for the Majorana neutrinos by employing a Breit-Wigner 
approximation, and taking the flavor matrix structure into account. This will 
allow us to identify effective pole masses and widths of the quasi-degenerate 
Majorana neutrino flavors. Then, we will use these propagators to determine 
a solution of the Kadanoff-Baym equation for the statistical propagator, which 
describes the non-equilibrium evolution of the Majorana neutrino fields. Finally, 
this solution will be used to compute the time-dependent lepton asymmetry 
produced by the interaction of the Majorana neutrinos with the thermal lepton 
and Higgs background. 



4.1 Retarded and advanced propagators 

The retarded and advanced Majorana neutrino two-point functions can be ob- 
tained by solving the algebraic Schwinger-Dyson equation (j66p in momentum 
space. Since we are interested in the time-dependence, we will also perform a 
Fourier transformation. 

In order to solve the Schwinger-Dyson equation (|66p. we insert the thermal 
lepton and Higgs propagators (|81l82p into the expression (|38p for the Majorana 
neutrino self-energy. The self-energy can be written in the form 

^N%A)iP) = -'2[ih^h),,PL + ih^hyPj,]Se^^^^^ip) , (85) 

where we have assumed even CP parities rji — 1 for simplicity. The function 
SiiPfi^j^^ip) denotes the lepton- Higgs loop integral. It can be written as a sum 
of vacuum and thermal contributions. Using dimensional regularization, the 
former is given by 

Si<^RiA)(p) = ^ (7 " + + 2 - In (^) ± *^e(p2)sgn(po)) p ■ 

The divergent contribution can be eliminated by introducing suitable mass and 
wave-function counterterms for the Majorana neutrino fields. If the difference 
between the Majorana mass parameters Mi^2 appearing in the Lagrangian is 
large compared to the decay width ri^2, it is convenient to use an on-shell 
renormalization scheme. However, in the present work, we are interested in 
mass splittings of the order or even smaller than the width. In this case, we 
find it more convenient to use MS renormalization. This means that the mass 
parameters are not necessarily equal to the physical pole masses Mf°2'^, and 
that the residues at the poles are not necessarily normalized. We will derive an 
equation for the pole masses in the following, and also properly take the residues 
into account in the computation of the lepton asymmetry. We stress that the 
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lepton asymmetry is independent of the renormalization scheme, and thus we 
can choose it to our convenience. 

Accordingly, we choose the wave-function and mass counterterms as 

= ^^(\~1E + H^^) + 2Y (87) 



1 
2 

In addition, it is convenient to choose the renormalization scale close to the mass 
scale of the Majorana neutrinos, such that the logarithm appearing in ([55)1 is 
small. We will set ^ = (Mi + M2)/2 for definiteness. We have checked that 
our numerical results are highly independent of the precise value in the range 
Ml < fi< M2. 

In a CF-symmetric configuration and for positive the thermal contribution 
to the lepton-Higgs loop is given by 



± 
i 



V 



± ^Q{p') / dIiidW^^{2^f5{p ^q-k)l^ [Uiqo) - Mko)] ■ (89) 



where, to shorten the notation, we have introduced p± = (±i?p,p) and q± — 
{±Eg,q). For negative p° the loop momentum k ~ (£'fc,k) in the second line 
of dHU) should be replaced with k — (E't, — k). For the analytical treatment, we 
will confine ourselves to the regime for which the Majorana neutrinos can be 
treated as non-relativistic when the lepton asymmetry is produced. In this case, 
the thermal contribution can be approximated by S£^*^f^j^^{p) ~ — r^|>/(12p^) ± 

2ip/{e^P°^^'^ — l)8(p^)sgn(po)/(327r). The first term describes a thermal mass 
shift, while the second one encodes the finite width. For the analytical treatment 
it is convenient to change the renormalization prescription in order to take the 
leading effect of the thermal mass shift into account, SZij — > SZij + SZij (T) 
with SZ,j{T) EE (/it/i),yrV(6/^^). The relation between the mass- and Yukawa 
coupling matrices at zero and finite temperature is 

MiT) = {PLZ{Tf + PRZ{Tf)M{T = Q){PlZ{T) + PrZ{T)*), 
ih1h){T) = Z{Tf{h^h){T^O)Z{Ty, (90) 

where Zij{T) = Vik{T){5kj + 6Zkj(T)). Here V{T) is a unitary matrix in flavor 
space that can be adjusted such that the mass matrix at finite temperature is 
diagonal and has real and positive entries, like the one at zero temperature. In 
the following, all quantities refer to the finite temperature values. We stress 
again that the lepton asymmetry is unaffected by the field rescaling, i.e. we 
would obtain identical results when setting SZij (T) to zero and using the vacuum 
parameters instead. The rescaling is performed for computational convenience. 
Altogether, when inserting the vacuum and thermal contributions to the self- 
energy as well as the counterterms into the Schwinger- Dyson equation (|66|. 
we obtain the following rcnormalizcd equations for the retarded and advanced 
propagators: 

[ {p - M,) <5*^- + ^(7,fePL + n^PR)p] S'^ (p) ^ -5'^ , (91) 
[{p Ah) S'' - *(7hPl + i:kPR)p] S'J ip) = ~6^^ . (92) 
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The coefficients are given by 
©(p^)sgn(|3o) 



7ij = ih^h)t. 



16tt 



+ i 



In^ 
167r2 



rp2 



rp2 

6^ 



(93) 

In the following, we will first identify the poles and the corresponding residua of 
the retarded and advanced propagators, and then Fourier transform the solution 
taking the contribution from the poles into account. 

Let us first consider the retarded propagator. In order to find a solution we 
make a suitable Lorcntz decomposition [501 EB HH] , 

Sr{p) = SllPl + SrrPr + SlrPl^ + SrlPr^ , 

where each coefficient is a two-by-two matrix in flavor space. We denote the 
expression in square brackets on the left-hand side of the Schwinger-Dyson equa- 
tion by r2, and perform an analogous decomposition: 



ik 
LR 

ik 



gik 



+ ilik , 



i2^^ = d + i^ki 



Then the solution of the SD equation reads: 

SrR = — [^RR — p'^^RL^j^^riLRl 



- p^Hlr^rr^rl] 



SlR — —^ll^LrSrr — IflRR^j^n^LL — 



RL 



(94a) 
(94b) 
(94c) 

(95a) 
(95b) 
(95c) 
(95d) 



Srl — —^rri^RlSlL — [^LL^Rj^^RR — P^^Lr] 

The components can be computed explicitly. It is convenient to write them as 

SXY 



S 



XY 



H 



(96) 



where X,Y E L, R, and using the common denominator H given by 

H = MiM2det[nRR-p^nRLnilnLR] 

= M1M2 det [^LL - p''nLRn-]^nRL] = Q V + M^M^ 

- p" [Af|(l -f i7ii)' + Ml{l + Z722)' - Ma/2(7?2 + ill)] 
Q = dot r^ifl = dot r^ijL = (I + «7ii)(l + 4722) +712721 ■ 
For the enumerators, we obtain 



sll 

SRR 
SLR 

Srl 



V 



(97) 



/ Mi(Af| - p\l + 1^22?) + Ahp^l^ 

ip2(Mi7i2(l + 4722) + M272i(l + i7ii)) 
zp2(7\/i7i2(l + ^722) + M272i(f + i7ii)) 

M^iMl -p\l + i7n)2) + Mip'^l^ 

5Ll|7i2^721 1 

M|(l + t7ii) -p2Q(i + 4-^22) »(MlAf27l2 +P^Q72l) 

i(MlM2721 +P^Qll2) M2(1 + ^722) - p2Q(l + i^ii) 

(98) 
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The denominator H can be written in the form 

H = Q^{p'' ~xi){p^ ~X2). (99) 

An the components are proportional to 1/H (note that H is invariant under 
the replacement 712 721)- Therefore, the retarded propagator has poles at 
~ a; 1.2. The two complex poles are given by 

{y±wf , 

^'■^ 4g2 Ip2=a;i 2 ' I 

where the quantities V and W have been defined as 



Y = yJ{M^{\ + 1112) - M2(l + i7ii))2 - MiM2{iv2 + l2if , 
W = \/{Mi{l + ^722) + Af2(l + ilii)Y - MiA/2(7i2 - 72i)2 . 

Note that these expressions are implicit equations for the poles, because the 
coefficients 7^- themselves depend on p^ . It is possible to determine the poles 
iteratively, starting with an initial guess that we choose as Jij{p'^ = /i^) with 
fj, = (Ml + M2)/2. The poles of the retarded propagator can be parameterized 
by effective dispersion relations ujpi and widths Fp/, 

p^ = ±uj.pi-'-Tpi, / = 1,2, (101) 

which are determined by = Pq — = x/, 

Vpi^ IM^/)I ^ ^p^^ /Rc(^,) + r2 /4 + p2 . (102) 

LOpI V 

The pole masses Mf"'"'^ = uipi \p=o and widths Pj"'*^ = Fp/jp^o obtained for 
vanishing spatial momentum are shown in Fig. [3] for a representative choice of 
parameters. We remind the reader that the masses Mi^2 are the renormalized 
mass parameters in the modified MS scheme where the mass shift due to the 
finite temperature is already absorbed into Mi 2- The relation to the vacuum 
mass parameters can be easily inferred from the rescaling prescription discussed 
above. 

It is instructive to derive approximate analytical expressions in two limiting cases 
of interest. Since we are mainly interested in mass splittings of the order of the 
width, (M2-Mi)/Mi ~ 0{{h^K)ij), and smah Yukawa couplings \{Wh)ij\ < 1, 
we will assume in all cases that (M2 — Mi)/Mi maxij |(/i'''/i),j |^. We also as- 
sume that the different entries are not fine-tuned, such that e.g. the difference 
{h)h)22 — {h^h)ii is of the same order of magnitude as the individual terms. 
Provided these assumptions are satisfied, we find that simple approximate ex- 
pressions can be obtained depending on the relative size of (M2 — Mi) /Mi and 
the real part of the flavor off-diagonal combination Re(/i^/i)i2/(87r) of Yukawa 
couphngs. Let us first consider the regime (M2 — Mi)/Mi > Re{h^ h)i2/ {^tt). 
Here we find that 



Mr" ^ M 



1 1 



F^ - r..i^M.(l + -,j^) . (103) 
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pole mass difference and effective widtfi 




(M2-AIi)/Mi 



Figure 3: Effective masses M^°^'^ = Wp/|p=o and widths r^"'*^ = ^pi\-p=o of the 
sterile Majorana neutrinos extracted from the complex poles of the resummed 
retarded and advanced propagators for {h)h)ii = 0.03, (/i^/i)22 = 0.045, 
{K^h)i2 = 0.03 • e*^/'' and T = 0.25A/i. The black dot-dashed lines show the 
approximation valid for (M2 — Mi)/Mi > Re{h^h) 12/(871), see Eqs. (|103|) . and 
the black dotted lines correspond to the approximate expressions (|104p valid for 
Re(/it/i)i2/(87r) > {M2 - Mi)/Mi > maji,j\{h'<h),j\^ . 



Thus, in this case, the pole masses are approximately equal to the renormalized 
mass parameters, and the widths arc given by the same expressions as in the 
hierarchical case. 

If, on the other hand, (A'h ~ Mi)/Mi < Rc(/itft.)i2/(87r), we find 

^ Ml + M2 ^ {M2-Ah){{h^h)22-{h^h)i,) 



2 2^i{h^h)22 - {h^h)^i? + 4(Re(/it/i)i2; 



2 



ole 



M, 
IGtt 



(l+ ^M./r_i ) {{h^h)^i + {h^h)22 



± ^J ((/it/i)22 - (/it/i)n)2 + 4(Re(;it/i)i2)2^ . (104) 



Note that the off-diagonal coupling parameter is bounded from above, | (/i^/i)i2 P < 
{h^h)ii{h'^h)22, such that the effective width cannot become negative. As can 
be seen in Fig. |31 the transition between the two regimes is relatively fast, i.e. 
the approximate expressions can be used for most of the considered parameter 
space. In addition, both cases agree if Rc{W h) 12 / {{h'^h) 22 — (^^^)ii) ^ 1 is 
small. We find that Eq. (jl03p yields a reasonable approximation within both 
regimes if the ratio is smaller than ^ 1/3. 

It turns out that the replacement 7^ — lijip^ = A*^) yields a rather accurate 
approximation for small Yukawa couplings and mass splittings of the order or 
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larger than the width. We will use this simplification in the following in order 
to obtain an analytical solution for the retarded propagator. In addition, we 
will assume that 1 + 2/(elP"l/"^ — 1) w 1 for simplicity. The resulting solution 
can be brought into the form 

Snip) - + 4^ , (105) 

where the residua Zir are matrices in flavor space, and can be decomposed 
analogously to the propagator, 

ZiR = ZjllPl + ZirrPr + ZiLRPLf + ZjrlPrP for / = 1, 2 . 

The components Zjxy {X, Y ^ L, R) of the residues can be computed explicitly 
using 

1 sxy{p^ =xi) 1 sxyjp^ X2) , . 

ZiXY ^ ' ^^XY = • (106) 

Q"^ Xl - X2 Q'^ X2 - Xl 

The solution for the advanced propagator can be formally obtained by replac- 
ing J7 — >■ (17)^, which implies VLxy ^xy- Using that ^ll and VLrr are 
hermitian, we find that the advanced propagator is given by 

Sa (p) - + sLPr + sIrPlP + Si^PRf , (107) 

where all expressions are evaluated at momentum p. For the poles, hermitian 
conjugation is equivalent to xi 1— >■ xj. Therefore, we obtain 

Sa{p) = ^ + ^ , (108) 
p'^ — xl — X2 

where Zja = z\j^j^Pl + 

z\llPr + zLrPlP + zIr^PrP for / - 1, 2. 
We are interested in the retarded and advanced propagators in time rather 
than frequency representation. Both representations are related by the Fourier 
transformation 

SR(A)^{t) - 1 ^ e-'^'«* Sr(a){p) ■ (109) 

The Fourier integration can be performed using the Cauchy theorem, and closing 
the contour in the upper imaginary plane for t < Q and in the lower one for t > 0. 
Note that, due to the non-analytic behaviour of at the threshold p^ — 0, it 
is necessary to use a contour which is obtained by dividing a semi-circle into 
three parts such that the lines corresponding to Rc(po) = ±|p| are left out. 
Within the analytical Breit-Wigner approach pursued here, we will take only 
the contributions from the poles into account, and neglect the contributions that 
arise from the parts of the contours parallel to the imaginary axis. This will 
be justified later on by comparison with numerical solutions. Since the poles of 
the retarded propagator have a negative imaginary part, the integral vanishes 
in the case t < 0. This is precisely what is expected for a retarded propagator. 
Analogously, the advanced propagator vanishes when t > 0. We find for the 
result 

Sr^A) pit) = ±&i±t) e^^-^'l-" [e-*--^*57^,(^) + e+^'^''^*^/H(^J (HO) 
7=1,2 
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where, the positive and negative frequency coefficients for retarded and advanced 
function, respectively, are given by 



± = ^l^lMjP^wj^^ (111) 



In the hmit of small Yukawa couplings and mass splitting larger than the width, 
the retarded and advanced propagators approach the well-known result 



S%A) p(i) ^ ±e(±t)eTr,,t/2 p7 gj^(^^^^) ^ .^^ cos(wp,t)j 5'^ . (112) 



The flavour off-diagonal contributions, which arc suppressed by the Yukawa 
coupling, are important for the creation of a lepton asymmetry. The leading 
contributions obtained by expanding in the Yukawa coupling are given in the 
Appendix. 

4.2 Lepton asymmetry 

The statistical propagator for the Majorana neutrinos, which describes the non- 
equilibrium time evolution due to the interaction with the thermal lepton and 
Higgs background, can be directly obtained by inserting the retarded and ad- 
vanced propagators into Eq. (|79p . Next, we obtain the lepton asymmetry gener- 
ated by the deviation of the Majorana neutrinos from equilibrium from Eq. (|80p . 
It is convenient to rewrite it in the form 

^dt) = /T^^trU^^' (0,0)<(i)-A4' (0,0)A'"=(t) 



(27r)3 

The lepton asymmetry depends on the initial conditions for the Majorana neu- 
trino, and on the time-dependent coefficients K which are given by 



Jo Jo 



In order to compute the asymmetry, we insert Eq. (|84p for the Icpton-Higgs loop 
Se,^p^{t" — t'). Using the result for the retarded and advanced propagators, we 
can now decompose the double-time integral. 



^ > ^ /,, 7=1,2 £fc=±l 

^p'^') ^ [(Sh^ ^ 7!7^,(K})Wt,K}), (113) 

^ 7,7=1,2 £fc=±l 



where all the are summed over + and — modes. The coefficients arc given by 

ifiii^k}) = -^{h^h),aoSf|PLS,i;;^PRStr^,, 

I'jAi^k}) = -*(/i^/i).,7o5!^:/PL5,47Pfi5^;^7o, (114) 
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and the time integrals are explicitly given by 



Lij{t,{ek}) 



dt' / dt"e 



where 



A = e2k + e-:>,q-eiUjpj + iTp,j/2, 



(115) 



(116) 



The double time integral can be easily performed now. Since the interesting 
dynamics occurs at time scales of t ^ l/Tpj, and the width of lepton and 
Higgs is much larger than the one of the Majorana neutrinos, it is reasonable 
to consider the limit where 

r^^ ■ t > 1 . (117) 



In this limit, 



Lij{t,{ek}) 



1 



2{A- B 
+ 



(A2 + ry4)(B2 + r,V4) 




(118) 



Lets identify which terms yield the leading contribution to the lepton asymme- 
try. Formally, we have to look for contributions which can be of order l/Fi, 
1 or 1/ (Ml — M2). Only the term proportional to \/{A — B) can be of that 
order, provided that ei = £4. Then for / = J wc get a contribution of order 
l/F/, and ioi I ^ J of order l/(ajpi — ujp2 +i{Tpi +rp2)/2). The latter will be 
important in the maximal resonant case. Thus, from now on we only take the 
1/{A — B) term into account. 

Furthermore, the terms involving A^ or B^ are unsuppressed only if the cor- 
responding real parts of A or B can become zero within the phase-space inte- 
gration. This is only the case for £2 = £3 = ei or £2 = £3 = £4, respectively. 
Altogether, this means the leading contribution arises when all £'s are identical. 



Llj{t,{ek}) L}j{t) ^ Lij{t,{ek})\ei=e2=e3=ei=e 

Then, wc find the following result: 



(119) 



LiAt) 



i^pi - t^pj + «e(rp/ + Tpj)/2 



Je{Upi-Upj)t-(Tpi+Tp,,)t/2 



(wpj -k-q + ieTp.,l2Y + TjjA (wp/ - k 



-^£Fp,/2)2 + F2 /4 



The entries on the diagonal can be written in a simpler form. They are real- 
valued and equal for both signs of £, L~j^j{t) = Ljj{t) = Ljj{t), with 



Li lit) 



1 - e^^p^* 



Rc 



F. 



{iupi -k-q + iTpi/2)^ + Tjj4 



(120) 



Note that the time-dependence precisely coincides with the time-evolution of 
the lepton asymmetry that would be expected from a Boltzmann approach. 



29 



The ofF-diagonal entries, i.e. for I ^ J, lead to an oscillatory behaviour and will 
be important in the maximal resonant regime. They fulfill the relations 

Note that all entries except Ln become suppressed in the hierarchical limit 
M2 ^ ujp2 ^ Ml. On the other hand, in the extreme degenerate limit M2 Mi, 
T2 — > Fi, all the entries Ljj become equal. 

Putting the results together, we obtain the following expression for the Icpton 
asymmetry: 



riLit) 



(27r)3 J (27r)32g^ 



(121) 



J=l,2 e=±l 



where 



and 



F' 



JI 



tr 



A5^'p(0,0)7S'=/-A5^'p(0'0m/ 



iff = IjMiek = e}) = -^{h^hyaoS'ZPLS,/;^PRS]^ho 



IjiiUk = e}) = ~^ih^h),,^oSfjPLSi/;PRSTno ■ (123) 



pp 



(122) 



Using the properties of L'jj{t), one can write the lepton asymmetry as a sum of 
Boltzmann-like and oscillatory contributions. 



TlLit) 



(pp f (Pq 



(27r)3 J (27r)32g 



FiiLii{t) + F22L22{t) 



+ F^2^e[LUt)) + Ff2lm(L2+ (0) 



(124) 



where the coefficients F contain all the information about the interaction, the 
CP violation and the initial state, 



F, 



II 



F, 



F^, 



II 



FTi , 

F12 + F21 



f: 



F^2 



12 — 12 ' 12 ^ 21 ^1 ^ 21 

t ( F-i o F-i n T -|- 



12 -^12 21 ^ 21y 
and the functions L{t) characterize the time-evolution. 



(125) 



Lii{t) 
Ltiit) 



1 



Re 



pi 



{ujpi-k-q + iTpi/2y+Tj^/4^ 



Tp + i{ui. 



•pi 



- tOp2) 

r£0/2 



(wpi-fc-g + irpi/2)2 + r2 /4 



(wp2 -k-q- iTp2/2y- + 



(126) 



where Tp = {Tpi +Tp2)/2. 

One non-trivial element that remains to be specified is the initial condition for 
the statistical neutrino propagator AS'fp(0, 0). In principle, we have a great 
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freedom to impose initial conditions, as far as the validity of the approach is 
concerned. Clearly, it is desirable to motivate the choice by some physical 
considerations. Concretely, we will assume that the neutrinos are in 'vacuum' 
at t = t' = 0, i.e. 

SFp{t,t%=t,^o = Sr;{t - t')\t=t'=o ■ (127) 

The Majorana neutrinos will then be produced by the Yukawa interactions with 
the thermal bath of Higgs boson and Icptons. This resembles a typical initial 
condition often used for leptogcncsis calculations. 

By definition, the above choice implies that the deviation from equilibrium at 
the initial time is given by 

ASFp{t,t%=t'=o = 5fp(< - t%=t'=a - S%{t - t%=t'=o ■ (128) 
4.3 Expansion in the Yukawa couplings 

By combining the results for the Icpton asymmetry (|124[) and the coefficients 
(|llip of the positive and negative frequency modes of the retarded and advanced 
propagators, it is possible to determine the resonantly generated asymmetry 
for a given set of Yukawa couplings {h)h)ij and mass parameters Mi^2- The 
approximations done so far allow to choose a mass difference of the order of the 
width, (M2 - Mi)/Mi ~ 0((/it/i)y )/(87r) < 1, which is the regime in which 
we expect a maximal enhancement. It is possible to simplify the analytical 
expressions considerably within the parameter region Re(/i^/i)i2 <C \{h^h)22 — 
(/it/i)ii| and (M2 - Mi)/Mi > Re(/i'f/i)i2/(87r). We wiU therefore first discuss 
the analytical solution in this regime, and then compare it to the full analytical 
solution in Brcit-Wigncr approximation as well as a full numerical solution. 
In the parameter region described above, the pole masses are approximately 
equal to the renormalized mass parameters Mi ^2 and the widths are equal to 
the usual expression for the decay rate ri^2 given in Eq. (|103p . The disper- 
sion relation and damping rate of the retarded and advanced propagators are 
given by ojpi ~ ^/Mf + p^ and Tpj ~ T jMi/ujpi. Then, it is possible to obtain 
explicit analytical expressions for the coefficients Fij defined in Eq. ()122p by 
systematically expanding the solutions for retarded and advanced propa- 
gators for small Yukawa couplings {K^h)ij <C 1. In particular, at leading order 
in the Yukawa couplings, the vacuum initial condition is given by 

A5i^ Ji,t')|t=t'=o ^ S^'IfdM ^^'^^ . (129) 

The result for the coefficients Fjj is shown in the Appendix, including also 
several intermediate steps of the calculation. The result obtained for the lepton 
asymmetry will be discussed in the following section. 

5 Result for the lepton asymmetry 

In this section, we will first compare our analytical result for the lepton asym- 
metry obtained within the Kadanoff-Baym approach in Breit-Wigner approxi- 
mation with the one obtained from a Boltzmann treatment. Then, we will com- 
pare various analytical, semi-analytical and numerical results obtained within 
the Kadanoff-Baym framework. 
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5.1 Comparison of KadanofF-Baym and Boltzmann 



Let us start with discussing the exphcit analytical result in Breit-Wigner ap- 
proximation, and compare it to the well-known Boltzmann approach. More 
specifically, we will first present the result that is obtained under the assump- 
tions described in section 14.31 By inserting the coefficients Fjj obtained at 
leading order in the Yukawa couplings into the expression (|124p for the lepton 
asymmetry, we obtain the following result for its time-evolution: 
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/i^D(Wp/)Rc 
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iMiTi + iM2T2 
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iMiTi + iM2r2 
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where q = |q| and k = |k| denote the Higgs and lepton energies, Upi = 
•yAMp + p^ are the energies of the Majorana neutrinos (/ = 1,2), and k ■ pi = 
kujpj—\q]. We have also defined Wp = 2ajpiajp2/('^pi+'^p2)- The time-dependence 
is described by the functions L{t). Neglecting Yukawa-suppressed contributions 
in the expressions (|126l) yields 



1 - e^^p^* 



[iopi -k-qf + ry4 



1 _ g-i(t^pi-<^p2)tg-rpt ^ 

k-q?+Tl^/^ 



Tp + i{wpi - ujp2 

r^0/2 

(wp2-fc-g)2 + ry4y ' 



where F^^ = Te + is the sum of the thermal width of lepton and Higgs, 
and Tpj = MjTj/uipi is related to the decay widths of the Majorana neutrinos. 
Furthermore, we have defined Fp = (Fpi -I- Fp2)/2. 

The flavor-diagonal contributions iii(i) and L22(i) feature a time-dependence 
that is also expected within the Boltzmann approach for the contribution from 
the neutrino species A'^i and N2, respectively. In addition, there exists a flavor 
off-diagonal contribution which has an oscillating behaviour with fre- 

quency given by the energy difference Wpi — Wp2 of the two Majorana neutrino 
species. The latter can be interpreted as the contribution from coherent transi- 
tions between the neutrino mass eigenstates A^i and N2 due to the off-diagonal 
Yukawa coupling Whi2. 

It is instructive to compare the result for the lepton asymmetry with the classical 
Boltzmann approach. Within the latter, the lepton asymmetry is determined 
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by the Boltzmann equation 

^ Boltzmann /• ^3 r ^3 r ^3^ , 



where 



dt {2'kYlOpi J (27r)32Q y (27r)32fc 

X ef^|A^jv,^,^|2 [/jv^(l + - /,) - (1 - 

1 lTa[{h'^h)]j] MiMj{Mj - Mf) 



{h^h)n 8tt (Mj - Mfy + {MjTj - MjTj)^ 



is the 'wave' contribution to the CF-violating parameter [3T] (see also [TU]). 
and \MNI-^^<|)\^ = 4fc ■ pi{h^h)jj is related to the tree-level matrix element for 
the decay Nj icfi. For the setup considered here, the time-evolution of the 
classical distribution function is simply given by fNi{t) = /fd(wp/)(1 — e"^"^*). 
In the following, we will compare the KadanofF-Baym with the Boltzmann re- 
sult in the hierarchical and in the degenerate limit of the Majorana neutrino 
spectrum. 

Hierarchical limit 

In the hierarchical limit M2 ^ Mi , the contributions proportional to the Fermi- 
Dirac distribution evaluated at the energy of N2, /fd(wp2), are exponentially 
suppressed for the relevant temperature range T Mi. In addition, the coher- 
ent contribution ij^(t) oc l/(rp + i{ujpi ~ Wp2)) ~ i/i^p2 is strongly suppressed 
compared to the flavor-diagonal contribution Lii(t) oc l/Fpi. The lepton asym- 
metry obtained in the Kadanoff-Baym approach is thus given by 



riLit) 



l-ai[{h^h)l2] Ml f d^p f d^q f d^k 



Sir Ah J (27r)3a;pi J (27r)32(7 J (27r)32fc 
X (27r)35(p - k - q) X ^ 



Ki-fc-9)2 + ry4 

1 - e^^pi* 



X 4fc -pi (1 -f /0(g) - fi{k)) fFoi^pi) 



For comparison, within the classical Boltzmann approach one finds that 
Im[(/it/i)22] Ml f d^p f d^q f d^k 



_, Boltzmann 



it) 



Stt M2 J (27r)3a;pi J (27r)32g J (27r)32fc 
X (27r)^(5(p — k — q) x 27r(5(wpi — k — q) 

1 — e^'^pi* 

X Ak-pi (1 4- f^(q) - ft{k)) /fd(wpi) 



Tpi 

Thus, the thermal width of lepton and Higgs F^^ = Ff + F^ leads to a re- 
placement of the on-shcU delta function in the Boltzmann equations by a Breit- 
Wigner curve, in accordance with [35] , 

2'K5{ujp — k — q) 



K-fc-g)2 + F2 /4 
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Degenerate limit 



In the limit where the masses arc quasi-degenerate, |Mi — M2I ^ Mi. 2, we ap- 
proximate Upi ~ ujp2 ~ ujp in aU terms except for those containing the difference 
of the energies. Then we obtain 



riLit) = 



MiAhiM^ ~ Ml) 



Stt {Ml - M2)2 + (A/iFi - M2r2)2 

d^p r (Pq r (Pk 

W^iJ (27r)32(? J (27r)32fc 



X (27r)'^(5(p - k - q) x 1 , 

K-fc-g)2 + ry4 

xAk ■ p{l + f^{q) - fl{k))fFD{<^p) 



(130) 



4 Re 



7=1.2 



p7 



2^ _ g-i(tJpi-wp2)tg-(rpi+rp2)i/2 

Tpl + rp2 + 2i(lLJpl - LUp2) 



For comparison, the Boltzmann result in the degenerate limit reads 



3oltzmann 



MiM2(A/| - Aff) 



Stt 



(M| 



d'^p 



MlY 

d^Q 



(MiFi - M2r2)2 

d^k 



(27r)3wpi y (27r)32(7 J (27r)32fc 
X (27r)^5(p - k - q) x 27r(5(wp - k - q) 



(131) 



x4fc • p (l-t- /0(g) - ft{k))fFD{ojp) 



E 

7=1,2 



1 — e^'^!'^* 



p7 



The first line of both results agrees, and corresponds to the resonant enhance- 
ment described by the 'usual' CP- violating parameter ef^. Note in particular 
that in the denominator the 'regulator' MiFi — M2r2 appears, in accordance 
with [21]. Furthermore, as in the hierarchical case, the on-shell delta function 
for the energies in the Boltzmann result is replaced by a Breit-Wigner curve in 
the Kadanoff-Baym result. 

Finally, the Kadanoff-Baym result features an additional contribution inside the 
square brackets in the last line as compared to the Boltzmann result. It results 
from the flavor off-diagonal contribution Lj^(t). As discussed above, it can be 
attributed to the contribution to the lepton asymmetry from coherent flavor 
transitions between the Majorana neutrino species. It becomes comparable in 
size with the Boltzmann-likc contributions when the energy difference Upi — ujp2 
is of the order of the average of the decay rates, (Fpi -I- Fj,2)/2, of the Majorana 
neutrinos. This case occurs when the mass difference AA/ = A/2 — Mi of the 
Majorana neutrino mass eigenstates is of the order of the widths F/. Such 
a situation has been frequently considered, because it leads to the maximal 
possible resonant enhancement of the CP asymmetry. 

The time-evolution of the lepton asymmetry nL(t,p), defined via nhit) = 
J d^p/ {2t:Y nLit, p), is shown in figure|l]for various choices of the mass-squared 
splitting A/| — Ml- When the splitting is very large compared to the product 
Af/F/ (/ = 1,2), the Kadanoff-Baym solution oscillates around the Boltzmann 
result, and the value at large times is very similar for both cases. When the 
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ni,(t,p)/nf°'*^™™"(i = oo,p) 




Figure 4: Time-evolution of the lepton asymmetry ni(t, p). Shown is the Boltz- 
mann result (dashed) and the result obtained from the KadanofF-Baym equa- 
tions, for various values of the mass-squared splitting M| — (solid lines). All 
curves are normalized such that the Boltzmann result approaches unity. For the 
widths wc have chosen Fi ~ O.OlA/i, F2 = 0.015A/i, and p = 0, for illustration. 



splitting is of comparable size as the product of the widths and the masses, there 
arc significant differences. Finally, when M| — < | A/2r2 — MiFi | , the lepton 
asymmetry obtained in the Kadanoff-Baym analysis is suppressed compared to 
the one predicted by the Boltzmann approach. The reason is that the coherent 
contribution has the tendency to cancel part of the Boltzmann-likc contributions 
to the lepton asymmetry. Note that to obtain figure S] we have neglected the 
effect of the thermal width of lepton and Higgs. Their inclusion would multiply 
the Kadanoff-Baym results for ni(i,p) by an additional overall factor, that is 
identical to the ratio of lepton asymmetries obtained from Kadanoff-Baym and 
Boltzmann analyses in the hierarchical limit. 

In order to quantify the suppression of the lepton asymmetry that results from 
the additional coherent contributions, it is useful to investigate the paramet- 
ric dependence on the mass-splitting, and to identify under which conditions 
the maximal possible resonant enhancement occurs. Within the Boltzmann 
approach, the resonant enhancement of the lepton asymmetry is maximized 
provided that the function 

r^Boltz^ann ^ MMMj - M^) 

- (Mi - MlY + (MiFi - A/2F2)2 ' ^'""^^ 

is maximally large. We may assume M2 > Mi without loss of generality. The 
maximum of the function Boltzmann occurs for a mass splitting of Af| — A/^ = 
AfiFi — A'jf2F2|, i.e. when the mass difference is of the order of the difference of 
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Figure 5: Function R characterizing the amount of resonant enhancement. 
Shown is the dependence on the mass-squared sphtting M| — for the Boltz- 
mann result (dashed) and for the result obtained based on Kadanoff-Baym equa- 
tions, for three fixed times t = 0.25/ri, l/FijOO (solid lines). The widths are 
chosen as in figure The dotted lines indicate the maximal resonant enhance- 
ment as predicted by the Boltzmann and Kadanoff-Baym results, respectively. 
The dot-dashed lines shows the numerical result obtained for large off-diagonal 
Yukawa couplings \ h''h\ 12 /{h^h)ii = 0.5 and 1, respectively. Here {h'^h)n are 
chosen as for the analytical result, and arg(/i^/i)i2 = 7r/8. The numerical and 
analytical results agree for |/i'''/i|i2 <^ \{h)h)22 — (^^/i)ii|, as expected. 



the decay widths of the neutrinos. The maximum is given by 

jDBoltzmann M1M2 /-. 

^max - 2|Miri -MaFal ' ^ ' 

Since the decay widths of the neutrinos are suppressed by the Yukawa cou- 
plings, this ratio can be much larger than one and can compensate the one-loop 
suppression of the CP- violating rates, so that it is possible to lower the scale 
of leptogencsis considerably below lO^GeV [TU]. In addition, it has been no- 
ticed that apparently a further enhancement occurs when not only the masses 
of the neutrinos arc quasi-degenerate, but also their decay widths Fi and F2 
are of similar size, so that Rmax « |MiFi — M2F2|~^ is further enhanced. Phe- 
nomenologically, this situation could occur when the degeneracy of the masses 
and of the couplings originates from an underlying flavor symmetry (see e.g. 

The Kadanoff-Baym result for the lepton asymmetry is proportional to the func- 
tion j^Boitzmann Howevcr, the additional coherent contributions to the 
lepton asymmetry have the tendency to cancel the Boltzmann-like contribution. 
Therefore, in order to determine the maximal possible enhancement, it is useful 
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to define a modified, time-dependent generalization by 



J^^B _ j^Boltzmann ^ 1—1.2 



J2 l~e~'^pJ* - 4Rc f l-e-'<"pi-"p2'*e-'^!'i+^f2)t/2 



rpi-i-rp2+2i(tijpi —ojp2) 



ik2 

(134) 

The dependence of i?^^ (i) and RBoitzmann mass-squared splitting - 

Ml is shown in figure [5l for a choice of parameters for which r2/Ti = 1.5. 
We observe that the resonant enhancement obtained in the Kadanoff-Baym 
analysis is smaller compared to the Boltzmann approach, in accordance with 
the above discussion. In order to identify the maximal possible enhancement, it 
is instructive to consider the value of R^^{t) for times t > l/Tpj, 

nKB(,\\ _ TjBoltzmann ^ (^pl - a^p2)^ + pi - ^ I 



j^Boltzmann ^ 



(M| - M2)2 + (A/iFi -I- M2r2)2 



(135) 



In the last line we have used that for quasi-degenerate masses Tpi = MjT i / copj ~ 
MjTi/Ldp and Wpi — ujp2 — {Mf — M|)/(2ajp). Thus, compared to the Boltz- 
mann treatment, the coherent contributions lead to an additional factor that 
effectively changes the regulator in the resonant term from the difference to the 
sum of the decay rates. This means the maximal enhancement of the lepton 
asymmetry occurs in fact for Af| — = AfiFi -I- M2T2, and is given by 

r>KB ^ M1M2 , , 

^marr 2(Afiri + Af2r2) ' ^ ' 

Consequently, the coherent contributions have the effect to cut off the resonant 
enhancement that occurs in the doubly degenerate limit Mi -> A/2 and Fi — r2 
within the Boltzmann approach. 

As is apparent from the expressions derived above, the maximal resonant en- 
hancement within the Kadanoff-Baym analysis is always smaller compared to 
the Boltzmann approach. The predictions of both approaches are comparable if 
either T2/ri ^ 1 or r2/Ti ^ 1, while the Kadanoff-Baym result is significantly 
smaller than the Boltzmann result when both decay widths are of comparable 
size. 



5.2 Comparison of analytical and numerical results 

The explicit analytical solutions obtained in the Kadanoff-Baym approach that 
have been discussed in this section have been obtained by performing a Brcit- 
Wigner approximation and an expansion for small Yukawa couplings that is 
valid within the parameter region Re(/i^/i)i2 ^ |(/i^^)22 — (^^^)ii| and (Af2 — 
A/i)/Afi >Re(/i'f/i)i2/(87r). 

We will now compare the analytical results with a semi-analytical and then 
with a numerical approach. The semi-analytical result is also based on the 
Breit-Wigner approximation. However, the coefficients Fij entering in the ex- 
pression (11241) for the lepton asymmetry are computed using the full solutions 
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for the effective masses and widths ujpj and Tpj given in Eq. (|102p as weh as 
the full expressions (^5)) for the retarded and advanced propagators. Conse- 
quently, this semi-analytical approach is valid also when Re(/i^ft.)i2 ~ {K^h)ii 
and for mass splittings (M2 — Mi)/Mi larger or smaller than Re(/i^/i)i2/(87r). 
We have checked that the results obtained in both approaches agree provided 
that \Wh\12 < Wh)22 - {,h^h)ii\. 

The semi-analytical result is shown as dot-dashed lines in Fig. [S] for two choices 
of |/i^/i|i2- The resonant enhancement becomes slightly stronger when |ft.^/i|i2 
is increased. However, the mass splitting for which the maximal enhancement 
occurs remains approximately the same. In addition, a good agreement is found 
also for extremely small mass splittings. Formally the validity of the analytical 
approach requires that (M2 - Mi)/Mi > Reih^i h)i2/ (Stt) - 0.005(0.01) for 
\h''h\i2/{h''h)ii = 0.5(1). Nevertheless, the comparison with the semi-analytical 
approach indicates that the asymptotic value of the resonant enhancement is 
predicted correctly also for extremely small mass splittings. It is worth noting 
that the semi-analytical result can be reproduced rather accurately by replacing 
the first line of Eq. (ITBUl) by 

Ini[(fet/^)2^] M,M2{Mi-Mf) 

Stt ((Aff'"')2-(A/f°'")2)2 + (Mf°'"rf'"'-M|'°'Tf'")2 ■ ^ ' 

and using the full expressions for ujpi and Fp/. Here the denominator contains 
the pole masses and widths discussed in section |4l and the enumerator is pro- 
portional to the basis-invariant quantity 

lm[tT[h''hM^MM^h^h*M]] = lni[{h^h)l2]MiM2{Ml - ) , 

related to CF- violation [5]. Since the pole masses and widths agree with the 
mass and width parameters in the domain of validity of the analytical result, see 
Eq. (jl03p . this replacement is consistent with the results discussed previously. 
The corresponding generalization of the enhancement factor leads to 



^l*>l/rp/ ((^Po'e)2 _ (Mf°'")2)2 + (^]^,jPolej.pole ^ j^jpole^ole-^ 
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This expression reproduces the numerical results shown in Fig. [5] with very good 
accuracy. It coincides with the expression obtained from the analytical solution 
for mass splittings that are large compared to Rc(/i^/i)i2/(87r), and also in the 
opposite limit, because then r^"'"^ -|- ppo''^ ~ Fi -f r2; as can be seen from 
Eq. (|104p . This behaviour is apparent also in Fig. \5\ 

In order to cross-check the Breit-Wigncr approximation, we have also solved 



the Kadanoff-Baym equation (j33bp for the spectral function of the Majorana 
neutrinos directly in time representation, for the spatial momentum mode p 0. 
The contribution from this mode to the lepton asymmetry can then be obtained 
via Egs. fSnirf^ . We have used the one- loop expression for the self-energy of 
the Majorana neutrinos, and inserted the thermal lepton and Higgs propagators 
(|81l82p . For the numerical calculation it is convenient to use an exponential 
momentum cutoff e^l'^l/^ for the computation of the Icpton-Higgs loop integral, 
and to solve the equation for the bare propagators. The renormalization is 
then performed by Fourier transforming the regulated loop integral to frequency 
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space, and applying a field renormalization such that the resulting expression 
matches with the renormalized Schwinger-Dyson Eq. (PT|) . For the cutoff we 
have used values in the range A/Afi ~ 3.3 — 5. We have checked that our 
renormalized prescription removes the dependence on the cutoff at a satisfactory 
level. Furthermore, we have used a time step at = 0.025/Mi and computed the 
solution for temperatures T ^ and T/Mi = 0.25 in the limit F^, F^ — > 0. For 
the renormalized mass and coupling parameters we have chosen the same values 
that are shown in Fig. 01 Wc find that the numerical solution for the spectral 
function can be well described by a sum of exponentially damped oscillating 
functions with damping rate and frequency in accordance with the Breit-Wigner 
solution (|110p . In addition, the resulting Icpton asymmetry agrees well with the 
analytical result for |/i^/i|i2 ^ |(^^ft-)22 — (^^^)ii| and with the semi-analytical 
result for Re(/i^/i)i2 ^ {h'^h)ii. 

6 Conclusions 

In this work we have studied the resonant enhancement of the lepton asymme- 
try that is generated by the lepton-number violating interaction of lepton and 
Higgs fields with a pair of quasi-degenerate Majorana neutrino species within 
a first-principle approach. The quantum-mechanical treatment of this system 
allows to study the resonant enhancement in the case when the mass difference 
of the heavy neutrino species is comparable to their width. This regime is of 
particular interest, because it corresponds to the parameter region for which the 
resonant enhancement is maximal, and thus determines the amount by which 
the efficiency of Icptogencsis can be increased compared to the case of a hierar- 
chical spectrum of the Majorana neutrinos. By treating the lepton and Higgs 
fields as a thermal bath and neglecting the backreaction, it was possible to ob- 
tain a formal analytical solution for the Kadanoff-Baym equations that describe 
the non-equilibrium evolution of the coupled system of Majorana neutrinos. We 
have used this solution in order to determine the generation of a lepton asym- 
metry, for which we have obtained approximate analytical as well as numerical 
results. Finally, we have compared the resonant enhancement that is obtained 
from the quantum-mechanical Kadanoff-Baym approach to the classical Boltz- 
mann approach. 

We find that for mass splittings of the Majorana neutrinos that are large com- 
pared to their widths the lepton asymmetries obtained within the Kadanoff- 
Baym and the Boltzmann approach have the same time-dependence, and differ 
by an overall factor that is related to the finite width of lepton and Higgs fields. 
When the mass splitting of the Majorana neutrinos is only slightly larger than 
their widths, the Kadanoff-Baym result is given by a sum of Boltzmann-like 
contributions and oscillating contributions which can be attributed to coherent 
transitions between the two Majorana neutrino species. If the mass splitting is 
of the same order as the width, a cancellation of the Boltzmann-like and coher- 
ent contributions occurs in the Kadanoff-Baym result and the lepton asymmetry 
is suppressed compared to the Boltzmann result. We have also compared the 
analytical approximation of the Kadanoff-Baym result to numerical results, and 
found good agreement within the domain of validity. In addition, wc find that 
the qualitative behaviour is robust. 

Our analytical result confirms the form of the regulator given in Eq. ([5]) . How- 
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ever, the cancellation that arises due to the coherent flavor transitions reduces 
the maximal possible enhancement. The analytical treatment indicates that 
the resonant enhancement can be estimated approximately by replacing the 
'regulator' according to MiFi — M2T2 -A/iFi + ]\I2T2- The cancellation is 
particularly important when the widths of the two neutrino species are of com- 
parable size. Wc therefore expect that the coherent flavor transitions arc crucial 
within phenomenological scenarios for which the masses and the decay rates are 
quasi-degenerate. However, since the interplay of Boltzmann-like and coherent 
contributions is a dynamical process, it is necessary to investigate these sce- 
narios in detail. Depending on the considered scenario, it may be appropriate 
to consider initial condition for the Majorana neutrinos that are different from 
the vacuum initial conditions considered here, which will also have an impact 
on the asymmetry. In addition, we note that it would be interesting to include 
the back-reaction and to implement the gauge interactions of lepton and Higgs 
fields in a systematic way. Nevertheless, we expect that the partial cancellation 
of the final asymmetry due to coherence effects is a generical result. 
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A Expansion in the Yukawa couplings 

In this Appendix approximate analytical expressions for the coefficients Fij 
defined in Eq. (jl22|) arc derived, under the assumptions discussed in section l473l 
These coefficients determine the resonant enhancement of the lepton asymmetry 
ni,(<)that is generated by Majorana neutrinos interacting with a thermal bath, 
according to Eq. (|124p . 

Within the parameter region described in section 14.31 the complex poles of 
the retarded and advanced Majorana neutrino two-point functions obtained at 
leading order in the Yukawa couplings arc given by Eq. (jl03p . 

XI = Mf-iMiTi, 7 = 1,2. (138) 
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For the residues of the poles we find in the same approximation, 



A/2 + ^ j + XI - X2 ^ ^ 2Ah [ 2A/2 + / 

.z(A/, + Af2)(A/2r,-A/,r2)-*^"('')^ « 



^2fl — 



4A/i V 

The positive and negative frequency coefficients for Fourier transformed retarded 
and advanced function can then be obtained using Eq. At zeroth order 

in the Yukawa couplings, the coefficients reduce to 

c,- ^ c,_ ^ i A/i + a;pi7o - p7 / 1 

1^ 2 Wpi ^0 

and the result for / = 2 is analogous. The leading contributions to the off- 
diagonal elements arise at first order in Wh and are given by 

12- _ A/i (A/i+Ljpi7o-p7)[Ar - ^Im(/it/ii2)(A/i - M2)75] 



St 



21- 



2ujpi M( - Mi - iMiTi + iM2T2 

A/i [AF - ji-lm{h^hi2)iMi ~ A/2)75](A/i + ujpgo ~ pj) 



il2 



^ 2u}pi Ml - A/| - iA/iFi + iA/2F2 

A^i (A/i - Wpi7o - P7)[AF - ^Im(/^t/^,2)(A/, - A/2)7.5] 



•-'IR — 


2wpi 




A/2 - A/2 + 


iA/iFi - iA/2F2 


q21 + 


Ml [ 


AF- 


- ^Im(/it/ii2)(A/i 


- M2)75](A/i - Wpi7o - P7) 




2wpi 
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iA/iFi - iA/2F2 


rl2- 


Ml 


(Mi 


+ Wpi7o - P7)[Ar - 


__i_Im(/it/i^2)(Mi- 1/2)75] 


'lA - 


2ajpi 




M2 - A/2 4 


- iA/iFi - iM2T2 
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[AF 


- ^Im(/it/ii2)(A/i 
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2LJpl 
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2wpi 
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- ^Im(/it/ii2)(A/i 
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2wpi 




A/2 - A/2 - 


- iA/iFi + iM2T2 
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Mz [Ar - ^Im(/itfei2)(Mi - M2)75](M2 + ujp2-fo - P7) 
2wp2 Ml - M| - iMiTi + iM2T2 

M2 {M2 + cjp27o - P7)[Ar - ^Im(fet/x,,)(Mi ~ 1/2)75] 
2wp2 - M| - iMiFi + iM2r2 

j\J2 [Ar - j^lTa{h^hi2){Mi ~ M2)75](M2 - u;p27o - P7) 
2wp2 Ml - Ml + iMiTi - iM2T2 

M2 {M2 - L0p2lo - P7)[Ar - ^Im(fet/^,3)(Ml - M2)75] 
2wp2 Ml - Ml + iMiTi - iM2T2 



S. 



12- 

2A 



S. 



2A 



^12+ 
^2A 



■'2A 



where 



M2 [Ar - ^lra{h^hi2){l\Ii - A/2)75](M2 + Wp27o - P7) 
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iMiTi - iM2T2 
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+ Wp27o - P7)[Ar - 


- jl^lm{h^h,2)iM, 


- ^2)75] 


2ujp2 
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- j^lMh^h^2){M, 


- ^2)75] (A/2 - Wp27o - P7) 


2iOp2 




Ml - Ml - 


iMiTi + iM2T2 




M2 


(A/2 


~ Wp27o - P7)[Ar - 


- ^Im(/^t/^l2)(A/l 


- A/2)75] 


2Wp2 




Ml - Ml - 


iMiTi + iM2T2 





Ar = (A/i + A/2) 



Rc(/i^/ii2) 



16n 



(141) 



The CP conjugated propagators differ by the sign of lm{h^ hi2). 

Using the vacuum initial condition (|129[) for the statistical propagator in leading 

order in the Yukawa couplings, the coefficients Fjj are given by 



Fjl = fFD{uJpl)tT 



A/i 



+/i=^r>(wp2)tr 



UJpl 

A/2 - 



'^p2 



ilji 



7,7/ ) 



(142) 



Lets consider for concreteness the first contribution. At leading order in Yukawas, 



7j/ 



-i^oAs%jPLSi/;^PRS],\^^o 



where 



AS'; 

AS' 



RI 



AI 



Re{h^h),2{Slli' - Sij) + ilm{h^h),2{Slli' 
Rc{h^h)i2{S\]'^ - S'^i) - tlm{h^h)i2{S\^' 



21e, 
Rl) ' 
12eN 
Al) ■ 



(143) 



(144) 
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From the explicit expressions for S^^^^^j given above, one obtains 

+ ^ . Im[(/itfe)^^] Ml {MiPl + M2Pb}[Mi - tOpila - P7] 
^1 * 167r ujpi Ml - M| + iMiTi - iMsFa 

Im[(/itfe)2^] Afi {MiPl + M2Pj^}[A/i + ^pi7o - P7] 
IGtt cjpi 1/2 - A/| - iA/iFi + iil/sFa 
. Ini[(fet/^)2^] M2 [M2 - cjp27o - P7]{AAPl + MaPfl} 

167r cjp2 Ml - M'^ + iA'/iFi - iM2r2 
^ lm[{h^h)\2] M2 [M2 + c^p27o - P7]{Mi^£ + M2PR} 
^2 * I67r Wp2 - M| - iMiFi + iM2r2 

. Im[(fet/t)2^] Ml [Ml - c^pi7o - P7]{MiPj^ + A-hPL} 

* 167r Wpi - M| - iMiFi + iA/2r2 
lm[{h^h)l2] Ml [Ml + c^pi7o - P7]{MiPfl + MzPl} 

* IGtt Wpi A/2 - A/| + iMiTi ~ iM2T2 
. Ini[(fet/t)2^] M2 {A/iPfl + A/2Pl}[A/2 - ^p27o ^ P7] 

^2 * iGtt ujp2 Af 2 - Af| - iAfiFi + iAf2r2 

lm[{h^h)l^] M2 {MiPr + M2Pl}[M2 + cjp2ln - P7] 



A^52 



ASX, 



AS 
AS 



A2 



167r Wp2 Ml - M| + iMiTi - iM2T2 



^11+ _ _ *r Afi -Wpi7o -P7 



2 ojpi 

oil- oil- «c Af 1 +^pi 70 -p7 
■5^^,/ = ^''i^ • (145) 

Z Wpi 



Finally, we obtain the following result for the coefficients Pjj: 

_ 7^- _ (fccjpi - kp)fi^{k,q)fFD{c^pi) Im[(fetfe)2^] 
^11 ^ ^11 



2fca;pi 47r 



A/JiM 

xRe 



^22 ^ ^22 



Ml ~ Ml + iMiFi - iAf2r2 ■ 

(fca;p2 - k.p)fe^{k,q)fFD{t^p2) Im[(fetfe)2^ 
2kujp2 47r 
M^ M2 

xRe- 



P± 



' Af2 - Af2 + iA//iFi - iA/2r2 ' 

(^pi +ojp2)fe4{k,q) Im[(fet/t)2^] 
4fcwpia;p2 47r 

1 J MiM2{kujp2 - kp)fFD{ujpi) MiM2{kujpi - \s.p)fFDiy!p2) 

2 1 Ml - A-f 2 ± iMiFi =F iA/2F2 Af f - Ml =f iAfiFi ± iM2F2 



X 



(Wpl + Wp2)/£0(fc, 9) Im[(/lt/l)22 
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X 



4fcti;pia;p2 47r 

1 J MiM2{ku!p2 - kp)fFD{ujpi) MiM2{kujpi - kp)fFD{iOp2) 

2 I Af 2 _ ^^2 -p jjVfiFi ± iAf2F2 Af 2 - Af 2 ± iAfiFi =F «Af2F2 
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